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T-H Using the AdS/CFT correspondence, we probe the scale-dependence of thermahzation in strongly 

coupled field theories following a quench, via calculations of two-point functions, Wilson loops and 
entanglement entropy in d = 2, 3, 4. In the saddlepoint approximation these probes are computed 
5_j in AdS space in terms of invariant geometric objects - geodesies, minimal surfaces and minimal 

volumes. Our calculations for two-dimensional field theories are analytical. In our strongly coupled 

^5 setting, all probes in all dimensions share certain universal features in their thermalization: (1) 

a slight delay in the onset of thermalization, (2) an apparent non-analyticity at the endpoint of 
thermalization, (3) top-down thermalization where the UV thermalizes first. For homogeneous 

T-H initial conditions the entanglement entropy thermalizes slowest, and sets a timescale for equilibration 

that saturates a causality bound over the range of scales studied. The growth rate of entanglement 

' ^ ' entropy density is nearly volume-independent for small volumes, but slows for larger volumes. 
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I. INTRODUCTION 



A. 



Motivation 



The observed nearly inviscid hydrodynamic expansion of the hot QCD matter produced in nuclear collisions at the 
Relativistic Heavy Ion Collider (RHIC) indicates that matter produced in these nuclear reactions is strongly coupled 
[iHS] . The inability of perturbation theory to account for this phenomenon has motivated studies of nonequilibrium 
dynamics in analytically tractable, strongly coupled gauge theories. The prototype of such theories is the maximally 
supersymmetric Yang-Mills theory in the large- iVc, large 't Hooft coupling limit, which is holographically dual to the 
classical limit of a superstring theory on the Anti-de Sitter space AdSsxS^ background Hydrodynamical evolution 
with a minimal shear viscosity emerges naturally as late-time behavior in the longitudinal expansion of a thermal 
gauge plasma in this model |S1 [5] . In the dual supergravity theory, the thermal state of the gauge theory is represented 
by a black brane in the asymptotic AdSs space, and the near-equilibrium dynamics of the gauge theory giving rise to 
hydrodynamic behavior is mapped onto the dynamics of perturbations of the AdSs-black brane metric. A nice recent 
review on the holographic study of hot QCD matter, focused on equilibrium and near-equilibrium aspects, is [7]. 

While the near-equilibrium dynamics of the strongly coupled super- Yang-Mills (SYM) theory is well studied, the 
process of thermalization itself is still poorly understood. The RHIC data demand that the time scale for equilibration 
of matter is considerably shorter than expected in the framework of perturbative approaches to thermalization |8l [9] . 
The dearth of other nonperturbative tools for the description of the short pre-hydrodynamic stage of the nuclear 
collision motivates the use of the AdS/CFT correspondence to study thermalization of strongly coupled plasmas. 
Phenomenologically, the central rapidity region of a high energy collision turns out to be boost invariant to a good 
approximation, so it is natural to consider boost invariant configurations [6l llOI - fT2] . One interesting question is how 
thermal equilibrium is reached in such systems. This is the question we will address in the present paper, albeit in 
the case of translationally invariant plasmas. 

We investigate the thermalization process not only in the strongly coupled (3-f 1) — 4— dimensional super- Yang-Mills 
theory, but also in the analogous lower (2 and 3) dimensional field theories. This has several motivations. The first 
reason is that the 2— dimensional version, which is dual to classical (super-)gravity on AdSa space, admits analytical 
solutions and thus allows us to explore thermalization over a wide range of parameters. A second motivation is that 
some version of the 2— dimensional dual theory, (1 -|- 1) dimensional super- Yang-Mills theory, can be formulated on 
a lattice and solved nonperturbatively jl3[ 114) . Finally, the comparison of theories in different dimensions makes it 
easier to distinguish generic aspects of thermalization of strongly coupled quantum field theories from aspects that 
are special features of 4— dimensional field theories. 



What is appropriately called the thermalization time may not only depend on the probe of the state of the field 
but also on the initial field configuration that evolves toward thermal equilibrium. Here one can distinguish between 
two broad classes of scenarios. The first class considers small perturbations around thermal equilibrium and asks how 
equilibrium is reached. For a gauge theory with a holographic dual, this approach amounts to studying the decay of 
small perturbations of an AdS black brane geometry. This is conveniently studied in terms of the quasi-normal modes 
of fields propagating in the black brane background [15]. The imaginary parts of the eigenvalues of these modes 
describe the thermal relaxation rates of various excitations, such as anisotropic perturbations of the stress-energy 



B. Holographic models of thermalization 



1. 



Thermalization scenarios 
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tensor. In the context of relativistic heavy-ion collisions, this approach is relevant for the study of viscous corrections 
to hydrodynamics and other transport processes in the presence of a thermal gauge theory plasma (see e.g. |16l for a 
review) . 

The second class of scenarios considers the thermalization of an initial field configuration, which is generally far 
from equilibrium. Through the AdS/CFT correspondence, the approach to thermal equilibrium in the boundary 
gauge theory is related to the process of black hole formation in the bulk. In early work, before the motivation from 
RHIC experiments, the goal was to understand how black hole formation from a gravitational collapse of a shell of 
matter would be encoded in the field theory in the boundary [17-20 . Also in more recent work, the initial conditions 
have usually been defined directly in the dual gravity theory - in most, if not all, of these cases the precise form of 
the initial condition is not known in terms of gauge theory excitations. Examples of such configurations are colliding 
gravitational shock waves |1T] , sheets of fundamental strings that may provide an AdS model of densely packed flux 
tubes [22] , and sudden perturbations of the metric near the boundary that propagate into the bulk [TI] [23] . Motivated 
by the properties of the initial state of a relativistic heavy-ion collision, which contains two highly energetic nuclei, 
these scenarios have in common that they inject energy into the AdS geometry at high momentum or short distance 
scales. 

Some of these recent studies of gravitational collapse in AdSs started from translation invariant, but locally 
anisotropic perturbations of the metric near the boundary and followed their propagation into the bulk, ultimately 
resulting in the formation of an event horizon [iF. The initial anisotropy of the metric is dissipated by the black 
brane, asymptotically resulting in viscous, boost-invariant scaling hydrodynamics on the boundary |121 124] . Gravi- 
tational collapse in AdS^+i induced by a scalar field perturbation was studied in [23], most explicitly for the case of 
AdS4, where the bulk equations of motion for gravity coupled to a massless scalar field were solved perturbatively 
for a small amplitude scalar perturbation that propagates in from the AdS boundary. The perturbation was taken 
to be translationally invariant along the spatial directions of the boundary and to vanish outside of a short time 
interval. An interesting technical result is that for this translationally invariant collapse in AdS4, the metric outside 
the infalling shell of matter coincides with a black brane metric to first non-trivial order in the perturbation [23]. 
This has the consequence that expectation values of local observables thermalize essentially instantaneously in the 
thermally quenched field theory dual to the infalling shell background. 

Similarly to the approach of 23 , in [25 the equilibrating field configuration was modeled by a homogeneous, 
infalling thin mass shell. The motivation for choosing such a configuration was that it naturally arises in the AdS 
dual description of the asymptotic limit of a transversely extended ensemble of fiux tubes |66| . The latter represent 
the energy density deposited by the colliding leading gauge charges (i.e. the valence quarks of the colliding nuclei), 
and the AdS description involves a sufficiently simple geometry to admit analytic or semi-analytic solutions. Our 
model does not allow us to study the approach to local isotropy and hydrodynamical behavior, but it enables us to 
investigate the approach of the field to a thermal configuration in momentum space and to measure the growth of its 
entanglement entropy, as discussed in Sections |II C| and |III C| In this sense, our present study is complementary to 
those mentioned above which focus on the approach to the hydrodynamical limit. Note that the work [25] considered a 
quasi-static approximation with metrics that correspond to stationary "snapshots" of the dynamic geometry containing 
the infalling mass shell. Here we study both the dynamic and the quasi-stationary geometries and show that there 
are some significant differences in the pictures that they yield of thermalization. Various aspects of the dynamic case 
of a thin infalling shell have been investigated in [26l - f28] . The effect of a dynamical boundary condition has also been 
studied in the background of a moving mirror j29| . New techniques for the evaluation of various Green functions in 
non-equilibrium backgrounds have been developed in 30], in addition to the complex contour approach by [31| . The 
counterpart to Chandrasekhar limit for a gravitational collapse of a degenerate star to a black hole in AdS/CFT has 
been investigated in |32j . Finally, thermalization after a quantum quench in gauge theory has also been modeled by 
a sudden change in couplings or other background fields than the bulk metric |331 134] . In probe brane approach such 
a change then can manifest as time-dependence in the induced metric on the brane. 

Concretely, we focus on the {d + l)-dimensional infalling shell geometry described in Poincare coordinates by the 
Vaidya metric 



where v labels ingoing null trajectories and we have set the AdS radius equal to 1. The boundary space-time is located 
at z = and x = (xi, . . . , Xd-i) correspond to the spatial coordinates on the boundary. For constant m{v) — M, the 
coordinate transformation 




(1) 



dv — dt 



dz 



(2) 



\- M z<^ 
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brings the metric ([ij in the form 



^ ' I- Mz^ 



dyi.'' 



On the boundary, the coordinates v and t coincide. We take the mass function of the infalhng shell to be 



m{v) 



M 
~2 



1 + tanh ■ 



(3) 



(4) 



where vq parametrizes the thickness of the shell falling along v — Q. We will often be interested in the zero thickness 
limit wq ~^ 0. The Vaidya metric describes a shell composed of tensionless null dust, which represents an analytical 
simplification of the tensionful shell models studied in _25 and of the model of [23]. 

The geometry outside the infalling shell (which corresponds to w > in the zero thickness limit) is identical to 
that of an AdS black brane geometry with Hawking temperature T = dM^/'^/Air, while the geometry inside the shell 
{v < when vq 0) is the same as that of pure AdS geometry. The causal structure of the Vaidya space-time is 
shown in Fig. [T] 



singularity 




AdS 
Schwarzschild 
{v > 0) 



FIG. 1: The causal structure of the Vaidya spacetime (in the vq — >■ limit) shown in the Poincare patch of AdS space. In this 
presentation, the asymptotic boundary (vertical line on the right hand side) is planar, and the null lines on the left hand side 
of the diagram represent the Poincare horizon. 



2. Probes of thermalization 

In order to probe the dynamics of the thermalization process, one can study a variety of observables of the boundary 
gauge theory. Expectation values of local gauge-invariant operators, such as the energy-momentum tensor and its 
derivatives provide valuable information about the applicability of viscous hydrodynamics, but they cannot be used 
to explore deviations from thermal equilibrium in detail. Nonlocal observables such as pair correlation functions, 
Wilson loop expectation values, and entanglement entropy provide much more information about progress towards 
thermalization. (They are also relevant to the physics probed in relativistic heavy ion collisions, e. g. through the jet 
quenching parameter q |35l I36j and the color screening length.) 

To illustrate, consider spatially homogeneous states of a weakly interacting massless scalar field. The energy- 
momentum tensor can be expressed as 

/(7k 
^fc^rn(k), (5) 

where n(k) = {a\M.-k) denotes the occupation number of a momentum mode of the scalar field. (There is no spatial 
dependence of the stress tensor because the field configuration is homogeneous here.) It is obvious from ([5| that 
this observable contains only limited information about the particle distribution. In particular, T^'' cannot inform 
us whether thermal equilibrium has been reached; it only tells us whether the pressure is locally isotropic. This is 
not even a sufficient condition for the matching of the field theory to hydrodynamics, because the equation of state 
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relies on the assumption of thermal occupation numbers for the various modes. It is also insufBcient to answer many 
other questions, e.g. about the expected spectrum of radiation emitted by collisions among the scalar particles. The 
equal-time two-point function, on the other hand, is given by 

— Kk) + l]exp(zk.x), (6) 

which allows us to extract detailed information about the particle distribution and thus probe its closeness to a thermal 
distribution in detail. Note that even the two-point function is insufficient to probe for the phase relationship between 
the occupation amplitudes of different field modes. This information can, however, be obtained by an analysis of the 
four-point function of the field. 

In an interacting quantum field theory, the two-point function is determined not only by the mode occupation 
numbers n(k) but also by the spectral density function a{k'^,'k). For the free field, (T(fc) ~ 2n5{k^ — m^), but more 
generally the spectral function is determined by the exact self-energy n(fc°,k) as follows: 

,(fcO ^ 2Imn(fc°,k) 

[(fc")2-k2 -to2 _Ren(fcO,k)]Vlmn(fcO,k)2' 

Since the self-energy is a function of the temperature, so is the spectral density function. For a weakly coupled field, 
the medium dependence is reflected, e.g. , in medium modifications of particle masses and widths. In the strongly 
coupled gauge theory, the vacuum spectral function has the generic form |37j 

CT(/fc°, k) - ((^0)2 - k2)''0((fc")2 _ k^), (8) 

where the exponent u depends on the dimension of the considered field. This becomes a sigmoidal function at finite 
temperature. The mode occupation number n(k) and the spectral function (T(fc°,k) combine multiplicatively to yield 
the true, physical excitation probability of a momentum mode k. The time evolution of the equal-time two-point 
function probes the approach of this excitation probability to thermal equilibrium. 

In the AdS/CFT correspondence there is a geometric intuition for why expectation values of local gauge- invariant 
operators are insensitive to details of the progress towards thermalization - they are only sensitive to phenomena near 
the AdS boundary. Thus they do not probe the details of phenomena occurring near the thermal scale. By contrast, 
nonlocal objects, such as two-point correlators of gauge-invariant operators and expectation values of Wilson loops, 
are dual to AdS quantities that probe deeper into the space-time and further away from the boundary. For example, 
at strong 't Hooft coupling, it is possible to approximate the path integral for the connected two-point correlator as a 
sum over all possible AdS geodesies connecting the two endpoints, which are placed on the AdS boundary [38]. (For 
certain subtleties regarding this statement, see [55WT] .') The geodesies probe the interior of AdS space, which is dual 
to the statement that they are sensitive to a wide range of energy scales in the boundary field theory. As we will 
see, Wilson loop expectation values 02 entanglement entropy [23 133 115 are also related to minimal lengths, 
surfaces and volumes of various kinds in AdS. These also extend into the bulk of AdS and hence probe a range of 
energy scales. 



C. Overview of this work 



We consider two-, three-, and four-dimensional field theories dual to gravity in asymptotically AdSa, AdS4, and AdSs 
spacetimes, respectively. The motivation for studying all these cases (rather than only the case of four-dimensional field 
theories potentially relevant for heavy-ion collisions) is that we are interested in generic conclusions on strongly coupled 
field theories with gravity duals. Lower dimensional theories are technically much simpler: analytic computations are 
possible in AdSs. Furthermore, our basic methods are similar in all dimensions - we calculate field theory observable 
that can be related to invariant geometric objects - minimal lengths, surfaces and volumes - in AdS spacetimes of 
different dimension. Our results show that many conclusions are not sensitive to the dimension. 

Sec. [n] discusses three measures of thermalization - two-point functions, Wilson loop expectation values, and 
entanglement entropy. We compute these observables at thermal equilibrium, and, in later sections, use deviations 
from these results as probes of incomplete thermalization. Sec. |III| relates thermal quenches in two-, three- and four- 
dimensional theories, to collapsing shells of matter in AdSa, AdS4 and AdSs that form black branes. The black brane 
formation represents the process of equilibration and thermalization of energy injected in the dual field theory. In 
these dynamical backgrounds we compute nonlocal observables as a function of time, and track the rate at which 
thermalization progresses for each observable at different spatial scales in the field theory. Sec. |IV| concludes the paper 
by summarizing our findings. An Appendix describes how the results would differ in a quasistatic approximation. A 
brief account of this work has appeared in [45j . 
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One broad finding is that these models exhibit a "top-down" thermahzation mechanism, see also [53]. We would 
like to stress that this is very much in contrast with the standard "bottom-up" paradigm [8] based on perturbative 
gauge theory. In the "bottom-up" scenario, hard quanta of the gauge field do not equilibrate directly by randomizing 
their momenta in two-body collisions, but they do so by radiating softer quanta, which gradually fill up the thermal 
phase space and equilibrate by collisions among themselves. In other words, infrared modes are the first to thermalize, 
and the thermahzation proceeds gradually ("from the bottom up") to more energetic modes of the gauge field. The 
bottom-up scenario is closely linked to the infrared divergence of the splitting functions of gauge bosons and fermions 
in the perturbative gauge theory. This contrasts with the "democratic" splitting properties of excitations in the 
strongly coupled SYM theory, which favor an approximately equal sharing of energy and momentum |461 14 7 j . We 
should therefore expect that the thermahzation process proceeds quite differently in the strongly coupled gauge theory. 
The phase space for radiation and, at strong coupling, also the spectral weight ([s]) is largest for the most energetic 
quanta. They then divide their energy and momentum approximately equally among their siblings, leading to a rapid 
cascading down to a thermal distribution. This consideration suggests that the thermahzation process more closely 
resembles a "top-down" scenario, in which energetic gauge field modes equilibrate first and soft modes last. Further 
broad conclusions appear in Sec. |IV[ 

It is appropriate here to comment on the differences and overlap between the current work and the recent related 
work in the literature. In |26l I27j. geodesies (or equivalently, entanglement entropy) in AdSa Vaidya were studied 
numerically for a thin but finite shell. In contrast, we also study geodesies in AdSa Vaidya analytically for an 
infinitesimally thin shell. In [28^ , Wilson loops (or equivalently, entanglement entropy) were studied in AdS4 Vaidya, 
for the case of circular and rectangular loops. We do not only confirm their results but also study geodesies in 
AdS4 Vaidya, which have not been closely investigated in the context of holographic thermahzation as far as we know. 
Moreover, we study AdSs Vaidya geometries for the first time by computing geodesies, circular and rectangular Wilson 
lines, and entanglement entropy for spherical volumes on the boundary. 



II. PROBES OF THERMALITY 



We are interested in processes where energy is injected into a strongly coupled field theory and equilibrates over 
time. Thus we need probes that take instantaneous snapshots of the state of the system which can be compared 
to the thermally equilibrated state. As discussed above, expectation values of local gauge-invariant operators are 
insufficient. Thus, in this section we describe three nonlocal probes - two point functions, Wilson loop expectation 
values, and entanglement entropy - which can be elegantly computed using geometric techniques in the AdS/CFT 
correspondence. We will show how these quantities computed in strongly coupled theory at thermal equilibrium 
differ from the corresponding quantities in the vacuum. In Sec. |III| we will use the same quantities as probes of 
thermahzation following a quench, by comparing the instantaneous values of the probes with the values at thermal 
equilibrium. 



A. Equal-time two point functions 

We begin our exploration with the equal-time two-point functions, which can be analytically calculated in two- 
dimensional conformal field theories at thermal equilibrium. We first remind the reader how the thermal occupation 
probabilities of momentum modes enter into the two-point function and how they are related to the spatial form of 
the equal-time Wightman function. We then explain how we generalize this approach to higher dimensional theories, 
where we can easily calculate the Wightman function in the saddle-point approximation. 



1. Two-dimensional free scalar field in thermal equilibrium 



We would like to understand how we can extract information about the thermahzation process from the time 
dependence of the Wightman function 



G^{t,x;t',x') = {0{t,x)0{t',x')), 



(9) 



where O is a local operator of dimension A. As a warniup exercise we consider the free massless scalar field in two 
dimensions with mode expansion 



(j){t, x) = 4>{x , x^) = 



dk+ 

k+>o 27r 



fe->0 



dk- 



(10) 
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We defined the operators as dimensionless, which imphes the commutation relation [afc+,a^,_|_] = k~^S{k^ — k'^), 
etc. Employing the free spectral function (Jo{k) ~ 27r(5(fc^) = 27r5(fc+fc~), this can be written as 



{2nf 

where for k'^jk^ > 0: 6^+ = o = '3^1+ 1 fc- = Qfc-j ^o,-fe- = ^I-- Consider the dimension one 

operators d~^(j) or d~(j). The vacuum two point function of is time and space translation invariant, so we consider 

G^(x±) = (a±0(O)a±0(x±)) = / ^(k^fao{k)e-^''^--^'-^^ = f ^fc+e-'^"-" ^ (12) 

For nonzero temperature, the spectral function of the free field remains unchanged, and we find 

G>(a;±;r)= / dk+ k+ (e-'^^''' {n{k+) + I) + e'''^''' n{k+)) , (13) 
Jk+>o ^ ' 

where n(/c) = {eP^ — 1)^^ is the standard (Bose) thermal occupation number. This can be rewritten as 

^k+g-^k+x+ ^+ (-^4) 

This is exactly the result expected from conformal invariance. The finite temperature computation is a computation 
on a Euclidean cylinder, which can be mapped to the complex plane using the exponential map, and we already know 
the answer on the plane. Undoing the coordinate transformation leads to the above two-point function. 



2. Strongly interacting scalar field theory in equilibrium from AdS/ CFT 



For an interacting field theory the spectral function changes with temperature. As a preparatory example of an 
interacting two-dimensional thermal field theory defined by a gauge-gravity duality we consider the AdSs case, i.e. 
the BTZ black hole, and make use of the results presented in |3S] . Writing w — (fc+ -I- fc^ ) /^/2 and k = (fc+ — k^)/ v^, 
the thermal Wightman function for a scalar field in such a theory is [48] : 



G>(fc) -z(l-coth— jim Gnik) 



-2i 



Im Gnik), 



(15) 



where Gu is the retarded Green's function and —2 Im Gr yields the spectral function. The precise form of the thermal 
retarded Green function depends on the dimension A of the field in the operator. For A = 1, we have [48] 



GniLj^k) 



1 

2^ 



1 



47rT 



1 



47rT 



Using Im ip{^ + iy) — ^ tanhTr?/, gives the spectral density of states as 



cr(w, k) — — 2Im Gr{u], k) 



sinh 



2T 



cosh ^ -I- cosh 



(16) 



(17) 



What can we learn from the equal-time Wightman function? One strategy is to use the vacuum subtracted 
Wightman function which at thermal equilibrium depends only on At = t' — t 



G>^^{At,k;T) = G>{At,k;T)-G>{At,k;T^O) = J ^ e"'"^* [G>(w,A:;T) - G>(c.,fc;T = 0)]. 



(18) 



Using the Bose occupation number density nsiuj) — (e"/-^ — 1) ^, the equal-time Wightman function {At = 0) is 
then 



G3>b(o,fc;r) 



= 2i 



did 

2^ 

= duj 
2^ 



nB(a;)o-(a;, k;T) + i 



duj 
2^ 
duj 
2^ 



d{-uj)a{uj,k;0) 



[a{uj,k;T)-a{uj,k-0)]. 



(19) 
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Inserting the explicit expressions, we obtain after some manipulations 



Gs>b(o,fc;T) 



2tt 



sinh 



2T 



cosh ^ + cosh ^ 



dw 



cosh ^ 



k 



2n cosh ^ + cosh ^ in 
Introducing the notation cosh a = coth(fc/2T), and using the integral identity 



did 

2T 



sinh a cosh h cosh a 

2T 





1 


+- cosh a 






2 


1 


— cosh a 



^ dLU 

2n 



k 

2T' 



(20) 



(21) 



the final result reads 



b (fc: 7^) = G> (0, fc; T) - G> (0, fc; T = 0) 



2tt e^lT-l 



(22) 



Having derived the functional form of the equal-time Wightman function for a thermally equilibrated state in AdSs, 
we ask what we could learn from the same function at non-equilibrium background. In that case time translation 
invariance is lost. As in the above, we could define a vacuum subtracted Wightman function G^,^(i, k) now depending 
on the equal-time t = t' , and expect that thermalization implies G^u\^(t,k) G^^^{k;T). This gives us a way to 
estimate equilibration times for different momentum scales k, by analyzing how this limiting behavior is approached. 
The result ( 22 1 , that the thermal contribution to the equal-time Wightman function of the two-dimensional field 
theory is proportional to the Bose distribution, justifies a more complete investigation of the equilibration properties 
of the Wightman function in two and more dimensions. For this purpose, we now turn to the methods made available 
by the holographic duality between strongly coupled supersymmetric gauge theories and gravity theories. We will 
also analyze thermalization times at different scales, in the spirit of the above motivation, but the details will differ 
somewhat. 



3. The geodesic approximation 

While the approach described above would give detailed information regarding scale-dependent thermalization, 
there is a technical challenge for strongly coupled field theories, even if they have a dual description in terms of 
gravity in an AdS space. The easiest case is two-dimensional field theory with a three-dimensional gravity dual. 
In this case, the thermal Wightman function in momentum space is known in closed form for operators of general 
conformal dimension A (A = 1 was analyzed above) |3H]- However, the integral with respect to uj that is needed 
to compute the equal- time Green function is sensitive to the ultraviolet completion of the theory and thus requires 
careful regulation. Thus, it would be more convenient to directly compute the equal-time thermal Wightman function 
from the AdSa theory. The next simplest case involves three-dimensional field theories with an AdS4 dual. In this case 
there is not even a closed form expression for the thermal Wightman function. Thus we cannot pursue the approach 
discussed above analytically. Furthermore, we are really interested in computing the equal-time Wightman function 
in a non-equilibrium setting in which matter collapsing in AdS to form a black hole models the thermalization of 
energy injected into the strongly interacting field theory. In this situation, there are no known analytic solutions for 
the Green functions in any dimension, and we have to resort to approximation schemes and/or numerical analysis. 

An insightful approach is to probe the non-equilibrium state of the strongly coupled field theory with a very heavy 
operator whose Wightman function can be approximated in terms of AdS geodesies as described below. According to 
the AdS/CFT dictionary, a massive scalar field ip{z,t,x.) with mass m in (d + 1) dimensions is dual to an operator 
0{t,x) of conformal dimension A{m,d) = | (d + + 4m2) [4 . In coordinates where the boundary of AdS is at 
z = 0, the scalar field ip behaves like 

(p{z, t, x) = z'''^cj,^o) (t, x) + . . . + (t, x) + . . . , (23) 

where the ellipsis corresponds to subleading terms in the z-expansion. In CFT language we interpret the non- 
normalizable mode (/>(o)(t,x) as the source to which 0{t,x.) couples and the normalizable mode 0(i)(t,x) as the 
vacuum expectation value of the renormalized operator 0^'°'^{t,x) 



0(i)(i,x)~(O-"(t,x)). 



(24) 
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The bare operator 0{t,x.) and the renormalized operator C'''''"(t,x) are related as [57] 

O(t,x) = zo^0-"(<,x), (25) 

where zq represents an IR bulk cut-off. The two-point function (©(t', x)C'(<, x')) is often computed using the on-shell 
(supergravity) action for if, when the action is explicitly given and solutions to the wave equation can be readily 
computed. 

However, for our purposes it is more convenient to follow |38j and to observe that the equal-time Green function 
can be computed via a path integral as 

(0(i,x)0(t,x')) =/"l?Pe*^^(^) « J2 e"^^- (26) 

geodesies 

The first expression sums over all paths that begin and end at the boundary points (t, x) and (t, x'), with L{V) being 
the proper length of the path {L(V) is imaginary for space-like trajectories). The second expression is a saddlepoint 
approximation to the path sum as a sum over geodesies ^S] (here C is the real length of the geodesic between the 
boundary points). The latter approximation is effective when the probe operator is heavy so that A ^ 1. It is easy, to 
check, for example, that in the zero temperature theory modeled by empty AdS space, this formula gives the correct 
conformally invariant equal-time two point function of the renormalized operator [38j 

(0^-(i,0)0'™(t,xo)) = ^. (27) 

Some care is necessary in evaluating Lorentzian correlators using the geodesic approximation, because careful consid- 
eration of steepest descent contours of integration is generally necessary [39ti4T| . 

The geodesic length C diverges due to contributions near the AdS boundary. Therefore, we define a renormalized 
length 5C = C + 21n(2;o/2), in terms of the cut-off Zq, by removing the divergent part of the geodesic length in pure 



AdS (see Eq. (51) below). The renormalized equal-time two-point function is 

(0(t,x)0(t,x')),.e„'-e-'^*^. (28) 

The geodesic approach to computing equal-time correlators gives a clear intuition for how thermalization proceeds 
in a strongly coupled theory with an AdS dual (Fig.[2|. Consider injecting energy homogeneously into the field theory 
- we model this in the AdS theory by dropping in a shell of matter from infinity. As the shell progresses inward 
into AdS space, the spacetime outside the shell will be well described by the AdS-black brane metric (except within 
the shell itself), while the metric inside the shell will be well described by the empty AdS metric. Now consider the 
two-point function calculation using geodesies. If the two points are close together on the boundary, the associated 
geodesic will not penetrate very much into the bulk space, and thus will not cross the matter shell. Thus, at these 
separations the correlation function should look thermal, because a thermal field theory is modeled by a black brane 
background in AdS. At larger boundary separations, the associated geodesic will penetrate the shell and be "refracted" 
by it (Fig. [2| leading to deviations from thermality. As time passes, the shell will penetrate deeper into the bulk and 
thus an ever larger range of spatial scales in the field theory will have associated geodesies that do not penetrate the 
shell and hence have thermal correlators. Thus we can come to a qualitative conclusion that thermalization proceeds 
top-down in this setup - ultraviolet, i.e. short distance, degrees of freedom thermalize first. In subsequent sections we 
will calculate the rates at which different spatial scales thermalize. 



4-. Two-point function for two-dimensional field theories in equilibrium: analytic computation 

In order to gain intuition for the two-point function in the geodesic approximation, we return to the case of two- 
dimensional boundary field theories. In order to obtain the thermal boundary-to-boundary Wightman function, we 
need to study space-like geodesies in the thermal black brane geometry in three dimensions. The metric is obtained 
by setting m{v) = M in ([T]), and is given by 

ds^ = -{r'-rj,)df + ^^+r'dx\ with t ^ v - iJ-^^^^ , r=-. (29) 

This is a black brane geometry with a horizon at r[{ = y/M |68| . The Vaidya geometry in three dimensions with a 
constant m{v) = M can be put in this form. 
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< ► 

Xo 



(A) (B) 

FIG. 2: (A) An example space- like geodesic that starts and ends on the boundary of AdS {z = 0) with a separation xo- Outside 
the shell, the geodesic propagates in a black brane geometry, while inside it propagates in an empty AdS geometry. The shell 
refracts the geodesic. The Wightman function at scales associated to geodesies that do not penetrate the shell of matter will 
be thermalized. (B) The minimal surface in AdS space associated to a circular Wilson loop. The shell of matter (indicated in 
light green) refracts the surface. Loops with associated surfaces that never penetrate the shell of matter will be thermalized. 
Both figures illustrate a quasistatic situation where the geodesic or minimal surface lies entirely at a fixed time. When the 
shell is dynamically falling into AdS, the geodesic or minimal surface, while remaining space-like, may not lie entirely within 
an equal-time surface. In both figures we are at late time when the shell is close to where the event horizon would be, so that 
the 'refraction' at the shell is clearly visible. 



Parametrizing the geodesic by the geodesic length A, the geodesic equations in the geometry (29) are 



1 = -{r'-r],)P 



, 2-2 



' H 



where E,J are constants and ' = d/dX. Plugging (301 and (31) into (32) and solving for r, 

1 



f = ±-\/r'* 
r 



(-1 + E^- J^yj^r^ + J^rjj 



Integrating this, we obtain 



or, equivalently, 



,A-An 



o e-(A-A„) 



(30) 
(31) 

(32) 



(33) 



(34) 



(35) 



Here, Aq is an integration constant which we will henceforth set to zero, Aq — 0, by absorbing it into the definition of 
the parameter A. Also, we defined the combinations 



A± = J2 - (1±£;)2, B± = {j±iy-E^. 



(36) 



The geodesic reaches the boundary r = oo (or z = 0) as A — > -l-oo. 

As we can see from ( 34 1 and ( 35 ) , the signs of A± , B± determine the behavior of the geodesic. The relation between 
the values of the parameters E, J and the signs of A±, B± is shown in Fig. [3j Note that not all possible combinations 
of signs occur. We discuss different cases (i)-(iv) shown in Fig. H in turn below. 

— 1/2 

(i) If A±,B± > 0, the geodesic is entirely outside the horizon. r(A) takes its minimmn value r = r^f[l + {A_^ + 
^i/2-j2^^ji/2 Xi-nin = (1/4) ln(j4+^_ ) while r(A) — t- oo as A ±00. One sample profile of this case is shown 

in Fig.gi). 
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FIG. 3: The values of E, J and the signs of A±, B±. The signs of A±, B± are written in the form (^4+ A-, B+ B-). See text 
for the behavior of geodesies in each regime of parameters. 



(ii) If A^A- < and B+B- < 0, the geodesic crosses the horizon r — th and reaches the singularity r = 0. Since 
in this case 



r(A) 



rH 



„2\ 



-A+A^ 



-2A 



2 



B+e-^) (p> 



> 



(37) 



(recall that we have set Aq =0), r(A) is a monotonically increasing function of A and the geodesic crosses the horizon 
only once. One sample plot of this behavior is presented in Fig. |4|^ii) . 
(iii) If A± < and B± > 0, the geodesic crosses the horizon twice at 



A 



1 



H± 



(38) 



but does not reach the singularity. One sample plot of this behavior is presented in Fig. Qiii). 

(iv) If A±,B± < 0, the geodesic crosses the horizon twice at (38) and furthermore hits the singularity twice at 



Ao± = -hi(-S±). 



(39) 



r is pure imaginary between Ao+ and Aq-- One sample plot of this behavior is presented in Fig.[4jiv). 

If one is interested in geodesies that start and end on the boundary of the BTZ black brane geometry ([29]), the 
only relevant cases are (i) and (iii). However, if one is interested in geodesies in the Vaidya geometry ([T]) for which 
the BTZ geometry (29) is only a part of the entire spacetime, then one should also consider cases (ii) and (iv) since 
they can appear as a part of geodesies in the full spacetime. 




-2 -1 1 

(i) £ = -0.1, J = 2 



2. 




1.5 




/0.5 




1. -0.5 


0.5 1. 


(ii) E = - 


-1.5, J = 2 




-3 -1 1 
(iii) E = -0.9, J = 0.01 



-0.5 0.5 1. 

(iv) £ = -5, J = 3 



FIG. 4: Sample profiles of r(A) (in bl ue) in the AdSs black brane background. The event horizon r_ff = 1 is shown in red. The 
integration constant Aq appearing in (34| has been set to Ao = — l/41n(|^+|/|j4_|). 



Substituting (34) into (30) and (31) and integrating the two equations, we find the expressions for t and x: 

.-A 



2rH 



oA 



=A 



A^e-^ 



(40) 
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where to, xq are constants of integration corresponding the values of x at A = +00. 



Comparing ( 40 1 with ([35]) , it can 
be seen that the quantity inside the absolute value of (40) vanishes or diverges at the horizon r = rn and thus t{X} 
diverges there. This is due to the well-known fact that the Schwarzschild time t is not well-defined at the horizon. 
Note that t{X) — > ^00 as A — > Xh± and therefore Xh+ {^h~) corresponds to the past (future) horizon. 
On the other hand, the Eddington-Finkelstein coordinate v defined by 



rHl 



2rH r[X) + rH 



(42) 



(see (29)) is well-defined across the future horizon A — Xh-, because the log divergence in t{X) gets canceled by the 



second term of (42) and ^(A) is finite across the future horizon. On the other hand, the divergence at the past horizon 
A = Xh+ is not canceled and v{X) diverges across the past horizon. The argument of the logarithm in (41), on the 
other hand, is always positive for r > 0; see (34). 



The above parametric expression of r, t, x, v in terms of A is useful for understanding the behavior of the geodesic, 
but for computational purposes it is useful to eliminate A and write t,x,v as functions of r. Inverting the relation 
(34), we obtain the following two branches for A(r): 



A±(r) - ^In 



1 



D{r) = + (-1 + E'^ ~ J^^h^ 



H- 



Note that we have been setting Aq = 0. Substituting (43) into (40)-(42) and setting = 0, we get 



t{r)-. 



2rH 



In 



x(r)± = — In 



v{r)- 



2rH 



In 



+ rH r^ + {E- l)rjj ± ^D{r) 



(43) 

(44) 
(45) 
(46) 



We will refer to the branch given by (i(r)+, x(r)-|_, w(r)+) as branch 1 and the one given by (i(r)_, a;(r)_, w(r)_) as 
branch 2. Only by combining both branches can we recover the full geodesic described by (l34|), (40) and (|4l|). 



geodesies end: 



From (40) and (41), we can compute the coordinate difference between the two boundary points on which the 



M = t{X -00) - t{X = 00) = 
£ = x{X = —00) — x{X = 00) = 



1 


log 


A+ 


1 


2^ 


A_ 


" 2^ 


1 


log 


B+ 


1 








2r^ 




" 2r^ 



log 



log 



il + E)^ 



j2-(i-£;)2' 

{J + lf-E^ 
{J-lf-E-^' 



(47) 
(48) 



We see that geodesies connecting endpoints at the same time (Ai = 0) correspond to taking E — {). If one is interested 
in such geodesies in the pure BTZ geometry, then we have only to consider the E = case. However, in section III 
we will be interested in geodesies in the Vaidya geometry, which are composed of two pieces stretching respectively in 
AdS and BTZ, and glued together at the shell location. In this case, one should also consider E ^ Q geodesies since 
they can appear as part of geodesies between equal-time endpoints in the full spacetime. 

Now let us focus on the equal-time geodesies with E = and compute the relation between the spatial boundary 
separation £ and the geodesic length £thcrinai in the black brane background. By setting E 
as 



in ( 48 ) , £ is computed 



rn J 



1 



1 



J — coth 



On the other hand, the geodesic length £thcrmai is computed from (43) as 



thermal 



A+(r = ro) - A_(r = tq) = 21n 



2ro 



^^W^\ 



2 In 



2ro sinh 



(49) 



(50) 



where = l/^o is the bulk IR cutoff and we dropped the ©(l/rg) quantity that vanishes as tq ^ 00. In the last 
equality, we used (49). 
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As discussed in Section II A 3 the geodesic length ( 50 ) is divergent in the tq — )■ oo hmit and should be regularized by 



subtracting the corresponding divergent part of the quantity in pure AdS geometry. Because the pure AdS geometry 
can be obtained by setting M — > or r// — > in the BTZ geometry (29), we can simply send r^/ — > in various 
quantities to obtain the corresponding quantities for pure AdS. By sending r^f — )■ in (50), we obtain the geodesic 
length in pure AdS: 



-Cpuro AdS = 21n(2ro) + (finite). 



Subtracting the divergent part of this from (50), we obtain 



sinh 

^'Cthormal(^) = 'Cthcrmal - 21n(2ro) = 2 In 



(51) 



(52) 



which allows us to compute the thermal renormalized two-point function defined in ( 28 ) . Observe that the two- 
point function result computed in the geodesic approximation coincides with the weak coupling result ( 14) for A = 1. 



Actually this is a consequence of conformal invariance, because the boundary CFT lives in noncompact space (namely, 
the X direction is not compactified in our setting). 



We can also obtain the expression for the equal-time geodesic in pure AdS geometry by sending rn — >■ in ( 45 ) . In 



doing this, we should send J — cx) at the same time because (49) becomes J = 2/{rHtj. In this limit, (43) and (45) 



give 



Ir 

A(r)± = ±cosh"^ — , 





(53) 
(54) 



where we dropped constants. The second equation (54) says that the equal-time geodesic is a semi-circle on the z-x 
plane, 

/IJ\ 2 

(55) 



Note that the pure AdS results (51), (53), (54) and (55) are correct not only for d — 2 but for general dimensions. 



because it is only the z-x plane that is relevant for such geodesies. 



Two 



-point function for d > 2 -dimensional field theories in equilibrium: numerical computation 



In any dimension we can consider a Vaidya type background ([T]) with constant mass function m(v) = M. This is a 
(d -|- l)-dimensional black brane geometry, and in the special case of three dimensions it can be put in the form (29 ). 
We consider space-like geodesies connecting the boundary points {t,xi) = {to,—£/2) and {t',x[) = {to,£/2) of such 
black brane geometries (for d = 3, we also have X2 = x'^ and for c? = 4, {x2,X'^ = {x'^^x'^). If we parametrize the 
geodesic in terms of the coordinate x\, which we will simply denote x in the rest of the Section, the solution will be 
given by the profiles v — v(x) and z — z{x). Inserting the IR bulk cut-off zq, the boundary conditions read 



z{-~ll2) = zo = z[ll2) , v{-ll2) = to - «(^/2) 
and are symmetric with respect to the z- and w-axis. The geodesic length reads 



thermal 



yj\ - 2z'{x)v'{x) - (1 - Mz{xY)v'{xY 

^1/2 



z{x) 



(56) 



(57) 



with ' = d/dx. The action is not explicitly dependent on x, which implies the existence of a conserved quantity. The 
conservation equation reads 

2z'v' - (1 - Mz'^y^ " (7) ^ ' 



1 



where is the value of z[x) at the midpoint. Notice also that w is a cyclic coordinate, implying a conserved 
momentum. From Eq. ([2|, it follows 

dz 

-^0-/ ^3^. (59) 
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Inserting the relation between v and z of Eq. ( 59 1 into the conservation equation ( 58 ) yields a first order differential 



equation in terms of z(x) only. We take the derivative of this equation with respect to x and solve the resulting second 
order differential equation 



(60) 



starting form the midpoint, where z(0) — and z'(0) = 0. We can then read off the boundary separation t 
corresponding to a specific midpoint value z■^, via z{i/2) — zq. By substituting the conservation equation into the 
expression for the geodesic length, we can simplify the expression to 



-'thermal 



dz 



^{\-Mz'^){l 



The renormalized length <5£thGrmai is 



(5 /^thermal (^) — -Cthormal (^) ' 



(61) 



(62) 



from which one obtains the renormalized equal-time two-point function through ( 28 1 . We plot the renormalized length 
'^'Cthcrmai as & function of the spatial scale £ in Fig. [5) 




FIG. 5: 5£thcrmai as a function of spatial scale I for d = 2 (red, dot dashed), d = 3 (green) and d = 4 (purple, dashed) for a 
black brane geometry with M = 1. The results for d = 2 agree with the analytical results of Section [II A 4| in the limit of a 
shell of zero thickness. 



B. Space- like Wilson loops 



Above we examined the Wightman functions for field theories at finite temperature in order to get a sense of 
what thermal equilibrium looks like according to this measure. A second non-local probe can be used to discuss 
the thermality of field theories is the (expectation value of the) Wilson loop. The Wilson loop is a gauge-invariant 
observable, constructed as the path-ordered contour integral over a closed loop C of the gauge field 

W{C) = j^TiiVefc^). (63) 

Wilson loops contain useful information about the non-perturbative behavior of (non-abelian) gauge theories, but are 
in general hard to calculate. In the AdS/CFT correspondence, the expectation value for the Wilson loop is related 
to the string partition function with a string worldsheet E extending in the bulk and ending on the loop C at the 
boundary 

{W{C)) ^ j VY,e~^^^\ (64) 

where we integrate over all inequivalent string surfaces E such that 9E = C at the AdS boundary and where A{Yj) 
corresponds to the string action. In the strongly coupled limit, we can consider a saddle point approximation for the 
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string partition function and reduce the calculation of the Wilson loop expectation value to determining the minimal 
area surface of the (classical) string worldsheet whose endpoints trace out the desired Wilson loop C on the AdS 
boundary 

{W{C)) ^ e-^^^^°\ (65) 

A{Tio) represents the area of the minimal area surface Sq with boundary C . The surface Eg is a solution to the 
equations of motion arising from the bosonic action of the string. Thus, the expectation value of a Wilson loop in 
the AdS/CFT correspondence is given by the area of a minimal surface in AdS space which is bounded by the loop 

c ga. 

In this section, we will compute space-like Wilson loop expectation values in the vacuum and in the thermally 



equilibrated theory (the black brane background) for d > 3-dimensional field theories. In Sec. Ill we will compute 
Wilson loops in a strongly coupled theory following a quench and compare them to the thermal results to assess 
progress towards thermalization. We will study two loop shapes - circles and strips as sketched in figure |6] - to test 
how and whether the shape of the loop affects the rate of thermalization. 



£ 




FIG. 6: Two Wilson loop shapes with their mmimal string surface: the circular Wilson loop (left) and the rectangular Wilson 
loop (right). 



1. Circular Wilson loops 

We first examine the circular space-like Wilson loop in a strongly coupled d > 3-dimensional field theory at finite 
temperature. At zero temperature the dual geometry is pure AdS^+i, and the associated minimal surface in AdS 
is simply a hemisphere |51j . This surprisingly simple solution to a complicated second-order non-linear equation of 
motion was obtained by taking an infinite straight Wilson line on the boundary, and applying a special conformal 
transformation to map this straight line to a circle. Since the conformal group of the boundary corresponds to 
reparameterizations of AdS we can apply the appropriate reparameterization to the AdS space, and it maps the 
original worldsheet - which was flat and extended straight into the AdS - to the hemisphere. This trick only works in 
pure AdS space where conformal invariance is preserved, so we are obliged to use numerical methods of solving the 
equations of motion when we want to calculate the Wilson loop for thermal or out-of-equilibrium states. 

For now, we want to study the minimal hemisphere surface in a thermal black brane background, with metric 
given by (IT]) with m{v) = M, where M is the constant tension of the black brane. On the boundary we choose a 
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two-dimensional plane {xi,X2) in which the circular Wilson loop is set. To parameterize the Wilson loop we introduce 
polar coordinates (p, f) . The minimal space-like surface with the circular Wilson loop as basis also has an azimuthal 
symmetry. The tip of the surface occurs at (v,z,x.) — (t;»,z»,0). The cross section at fixed z and w is a circle, and 
thus the surface is conveniently parameterized in terms of the radii p of these circles. The Nambu-Goto action for the 
string with circular symmetry in the {d + l)-dimensional AdS black brane background is 

^NG = ^ r dp^Jl- (1 - Mz{pY)v'{pY - 2z'{p)v'{p), (66) 

where we integrated out the (ys-factor as both z and v are independent of (p due to the circular symmetry. This action 
has an explicit p-dependence, implying that the second order differential equation in z{p) cannot be integrated to 
a first order differential equation. The coordinate v is still a cyclic coordinate, which implies that the associated 
momentum is conserved. To solve for an equal-time Wilson loop, we set the conserved energy to zero, which implies 



a relation between v and z as given in cq. (59 1. Inserting this relation in the action. 



the resulting equation of motion reduces to a second order differential equation in z{p) 

2 z' 1 \ ,^ . , ^ 12 



(67) 



z ' pi- Mz'i J^^ Mz +z ) ^ ^ _ . (68) 

One can easily check that the hemisphere is an analytic solution in pure AdS (obtained by setting M = in the 
above). It is given by 



z{p) = y^F^. (69) 

In case of a black brane background (M 7^ 0) we have to resort to numerical means. We impose boundary conditions 
at the midpoint assuming reflection symmetry along the z-axis 

z(0) = z,, z'(0) = 0, (70) 



and read off the boundary radius z{R) = Zq. Note that because the second term in the first bracket of (68) causes 



numerical issues at p = (the numerator and denominator both go to zero, but the ratio should take a fixed value). 



we cannot construct the solutions starting from the midpoint. Therefore we choose to solve (68) in the neighborhood 



of the midpoint by expanding around p — to quadratic order (odd powers vanish by the symmetry) 

Zpip) = — p . (71) 

The boundary conditions for the numerical solution are obtained by matching at a point pp close to the midpoint; 
z{pp) = Zp{pp) and z'{pp) = z'p{pp). 

The logarithm of the expectation value of Wilson loop operators is approximated by the area (-Athormai = a'^Nc) 
of the minimal surface. The largest contributions to the area are coming from the near-boundary region, because of 
the diverging AdS volume there. We subtract the divergent piece, which is proportional to the UV-cutoff I/zq. This 
procedure is interpreted in the dual CFT as removing the UV divergence of the Coulomb self-energy of a point charge. 
Subtraction of R/zq gives the regularized area 

R 

^-4thcrmal(^) = -^thermal (-R) • (72) 

The regularized area is depicted for AdS4 and AdS 5 in Fig. [7] 



2. Infinite Rectangular Strips 



A less symmetric Wilson loop is the rectangular strip parametrized by the boundary coordinates {xi,X2)- Assume 
that the strip is translationally invariant along the a;2-axis, such that the profile of the associated minimal surface 
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^■^thermal 




FIG. 7: 5.4thcrmai against the boundary radius R. The curves represent the string surface areas in a black brane background 
with unit mass for d = 3 (green) and d = 4 (purple, dashed). 



in AdSd-|_i is described by z(xi) and v{xi). In the following, we again denote x = xi. We impose the boundary 
conditions 

z{~£/2)^zo^z{£/2), v{-£/2)^to^v{£/2). (73) 

As for the circular Wilson loop, the equation of motion for the minimal surface is obtained by minimizing the Nambu- 
Goto action. In this case, because of the symmetries, the analysis closely resembles that of geodesies with some small 
differences in the action and the equations of motion. The area for a segment stretching between X2 G (0, R) in 
Eddington-Finkelstein coordinates is given by 

R f'/\ yi - (1 - Mz{xr)v'{xy ~ 2z'{x)^) 

The action does not exhibit an explicit x- or u-dependence, hence we can use the existence of two conserved 
quantities to simplify: these can be taken to be the midpoint of z(x), z*, and the momentum conjugate to v{x). The 
z„ conservation equation reads 

\~{\-Mz'')v'^ ~2z'v' =(-\\ (75) 



while the conservation of the momentum conjugate to v{x) implies the same relation as in Eq. (59). Inserting this 
relation in the conservation equation gives a first-order differential equation in z(x). We then take the derivative with 
respect to x to obtain a second-order differential equation in z(x) 

z" ^-2'^ + ^—^Mz''~^zt + -Mz''-\ (76) 

We impose the boundary conditions z' [Q) = 0, z(0) — at the midpoint, and read off the boundary separation i via 
z{£/2) = zq. Inserting the conservation equation in the Nambu-Goto action we get the following formula for the area 
(-^thermal = a'^Nc) of the minimal surface 

-Athormal — — "J" — , = • ( ' ' ) 



TT 



1 - Mz" 



The regularized area reads 



1 R 

<5Ahcrmal(^, R) — -^thermal , (78) 

Zq TT 



and is plotted in Fig. [8] for a four- and five-dimensional background. We have chosen to plot the results after dividing 
out R/n as we are concerned purely with the dependence of the area on the boundary separation £. 
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FIG. 8: 5.4thormai/(JJ/7r) as function of the spatial scale I. The curves represent the string surface areas in a black brane 
background with unit mass. The base of the surface is a rectangular strip at the AdS boundary. The green curve corresponds 
to d=3 and the purple, dashed curve to d=4. 



C. Entanglement entropy 

It is still an open question whether a suitable notion of "local entropy" density exists which is valid out of equilibrium 
and which satisfies some basics physical properties; in particular it must be a non-decreasing function of time. To define 
a local entropy density using gauge/gravity duality, one typically employs horizons since their area increases in time, 
and it has been suggested in 11 , 52^ that apparent horizons provide more compelling notions of local entropy density 
then global horizons do. In our case, the entropy associated to the apparent horizon instantaneously thermalizes, 
in the limit of a zero thickness shell, and does not provide a useful probe of thermalization. Skepticism against the 
identification of entropy with apparent horizon area in full generality has also been raised in a number of other papers, 
including HZlElliS. 

Another way of assessing thermalization at different spatial scales in an out-of-equilibrium theory is by measuring 
the entanglement entropy associated with volumes of different sizes and shapes. To review, consider dividing a 
quantum mechanical system in a state |$) into two spatially disjoint parts A and B. The density matrix of the 
system restricted to A and its entropy are computed by taking the traces 

PA -trB|$)($|; S{A) = -\.xa{pa\^Pa)- (79) 

In a (1 + 1) -dimensional conformal field theory, where A is an interval of length the entanglement entropy can be 
calculated analytically, yielding the universal result [SU |SS] : 

c ( l\ ... c f (i . Tr£\ 

vacuum : S'o(^) = - In — ; thermal equilibrium : St{(-) = - In — sinh — . (80) 
6 \a J 3 yna p J 

Here a is the UV cutoff of the field theory, c is the central charge, and /3 — T^^ denotes the inverse temperature. 

The entanglement entropy S{A) describes the amount of information loss associated with the restriction of an 
observer to the volume A. In the vacuum state for d > 2 dimensions, So{A) is proportional to the surface area of 



A [56j[57]; in d = 2 dimensions, as Eq. (80) shows, Sq depends logarithmically on the length of the interval £. At 
non-zero temperature, S{A) receives an additional contribution, which can be interpreted as thermal entropy [43l I44j . 
In the limit T — >■ cxd, the thermal contribution is proportional to the volume of the region A, just like the statistically 
defined thermal equilibrium entropy. Computing the time dependent entanglement entropy as a function of spatial 
scale and studying its approach to St{£) thus will provide a probe of scale-dependent thermalization. 

There is a precise proposal for computing entanglement entropy in strongly coupled field theories with AdS duals 
PSI HI] , where we think of the field theory as living on the boundary of an AdS space. Specifically, consider 
the boundary dA of a connected region in the field theory whose entanglement entropy we wish to compute. For a 
two-dimensional field theory, OA is a pair of points, for a three-dimensional field theory OA is a closed curve, and for 
a four-dimensional theory a surface. Now construct the minimal surface aA in AdS space that meets dA on the AdS 
boundary. For a two-dimensional field theory, aA is a geodesic in AdSs that approaches the boundary points dA, 
while for a three-dimensional field theory, a a is a minimal surface in AdS4 with boundary dA. For a four-dimensional 
field theory, a a is a minimal volume in AdSs which ends on the surface dA on the AdS boundary. The entanglement 
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entropy of the region A of the field theory is then given by 

S. = (81) 

where Gat denotes Newton's gravitational constant (and "Area" stands for the length of geodesies, the area of 2- 
surfaces, and the volume of 3-surfaces). The authors of [26l|43l|44] showed that this formula precisely reproduces the 
universal entropy formula of two-dimensional conformal field theories in thermal equilibrium |54ll55j . The restriction of 
the entanglement entropy to a finite spatial volume can be understood as a kind of coarse graining, by the discretization 
of the available momentum space modes, of the information contained in the quantum state |58j . 

Thus, to measure entanglement entropy in two dimensional field theories, we need to calculate space-like geodesies 
in the dual asymptotically AdSs background. This means there is an intimate relation between two-dimensional 
entanglement entropy and the equal-time Wightman functions computed in the geodesic approximation. Indeed, the 
two quantities can be interpreted in terms of each other, as was also observed in |27j . For three dimensional field 
theories, we need to compute space-like, two-dimensional minimal surfaces in a dynamical, asymptotically AdS4 back- 
ground. These coincide precisely with classical space-like worldsheets of fundamental strings with a fixed boundary, 
as previously studied in [5S]. As we discussed above, the exponential of the action of such objects (essentially the 
area) is known to give the expectation value of Wilson loops in the field theory |42]. Thus, in these three-dimensional 
field theories, there is an intimate relation between entanglement entropy and Wilson loop expectation values, as we 
explore in this article. 

For four dimensional field theories, the proposed measure of entanglement entropy for spherical regions involves 
minimal volume hyper-surfaces whose cross-sections with respect to the AdS radial coordinate are spheres. In par- 
ticular, in AdSs we can look for a 3-dimensional volume with 2-sphere cross-section which terminates on the AdS 
boundary. This allows us to examine the thermalization of the entanglement entropy in the 4-dimensional field theory. 
The methods used to compute these volumes are very similar to those for the circular Wilson loop described in Section 
|IIB H with only some powers changing with the dimension. Briefly, the action for the general spherically symmetric 
p- volume in AdS^ [p l£ d — 1) black brane background is 

Vthcrmai = As.-^ dp^^l - (1 - Mz{pY)v'{pY - 2z'{p)v'{p), (82) 
where Agv-i is surface area of a unit 5*^^^. Things continue to proceed as for the Wilson loops: the solution in the 



pure AdS {M = 0) background is now a hemi-hypersphere and still obeys ( 69 1 . What changes is the volume evaluated 
on the solution, which we evaluate with a cut-off near the boundary and subtract to renormalize as before. The result 
in the vacuum {M = 0) is 

VAdS = 4^i? / dpj-^'^^ = 2^ ^ + log + (finite) . (83) 







The regularized volume is calculated analogously to previous cases by subtracting the divergent part of the volume 
of a pure AdS 3-surface with terminates on the same sphere on the boundary, 

^Vthcrmal(i?) - Vthormal - 2^ f ^ + log . (84) 

We plot the regularised volume (5Vthcrmai as a function of the radius R of the sphere on the boundary in Fig. [9j 



III. SCALE DEPENDENT THERMALIZATION FROM ADS/CFT 

In the previous section we discussed how the equal-time Wightman function (related to geodesies in AdS space) 
and Wilson loops (related to minimal surfaces in AdS space), and entanglement entropy (related to minimal lengths, 
surfaces and volumes in different dimensions) can probe thermal equilibrium. In this section we will use the same 
probes to study the dynamics of equilibration in strongly coupled field theories. The basic setup is to drop a shell of 
matter with vanishing rest mass ("null dust") into AdS space to form a black hole. This is dual to the homogeneous 
injection of energy into the dual field theory and its subsequent thermalization. We want to probe the rates at which 
thermalization occurs on different spatial scales. To this end we consider dynamical Vaidya-type backgrounds, again 
given by equation ([T]), but this time using the continuous mass function of equation Q. This dynamical background 
continuously interpolates between pure AdS and AdS with a Schwarzschild black brane. Recall that the parameter 
vq in the mass-function Q dictates how "thick" the shell is. f g — is the step function limit corresponding to an 
infinitely thin shell. In this limit, the infalling shell represents a shock wave. 
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FIG. 9: Regularized volume iJVthermai as a function of the spatial scale R for d=4. 



A. Equilibration of the two-point function 

1. Two-dimensional field theories: analytic treatment 



We first consider the d = 2 case where analytic computations are possible. In this case, the Vaidya metric ([T]) 
becomes 



[r^ — m{v)]dv'^ + 2drdv + r^dx^ , r ; 



(85) 



In order to have analytic control over the solution, let us consider the thin shell limit — > of the mass profile Q, 
which gives 



m{v) = r'jj9{v) 



TH = VAf, 



(86) 



where 9{v) is the step function. 



Outside the shock wave, t; > 0, the metric is the standard planar black brane metric (291 

dr^ 



dsi 



-{r' - rif)dt' + 



+ r'^dx'^ 



t 



• H 



^in^:^ 

2rH r + rn 



(87) 



which we studied in detail in the previous section. Inside the shock wave, v < 0, the metric is the Poincare AdSa 



ds,„ = 



-r^dt^ 



dfr^ 



r^dx^ 



t = v- 



We would like to study geodesies in the AdSa Vaidya geometry which is (87 1 and (88) glued together across the 
infalling shell at w = 0. In particular, we focus on the equal-time geodesic which starts and ends at the same time 
i = ^0 on the boundary r — oo. There are two possible cases: (i) the geodesic does not reach the shell and is entirely 
outside of it, and (ii) the geodesic crosses the shell. Because case (ii) is more involved than case (i), let us first discuss 
case (i) briefly and then turn to the discussion of case (ii). 

In case (i), the geodesic is given by the equal-time geodesic in the pure black brane geometry, namely by (44|-(46) 
with E = 0. Therefore, the relation between the (renormalized) geodesic length i5£thermai a-nd the boundary separation 
i is given by ([52|. 

Now let us turn to case (ii). In this case, the part of the geodesic that is inside the shell is given by the equal-time 
geodesic in the pure AdS geometry, namely by (54), which we write in the following form; 



x{r) 



(89) 



where we set = 2/£. On the other hand, the part of the geodesic that is outside the shell is given by the geodesic 
in the pure black brane geometry, namely by (44)-(46). In the present case, we should not set the parameter E to 
zero, because the geodesic gets refracted at the shell and it does not have to be in a constant t slice. 



21 



The two parts of the geodesies, the inside part and the outside part, should be connected so that the total geodesic 
length is minimized. Just like Snell's law, this can be stated as a refraction law for the angles entering and exiting 
the shell, as follows. Let us write the metric (87 1, (88) as 



ds^ = -f(rfdv^ + 2drdv + r'^dx^ 



where 



H 



V <0 
w > ' 



(90) 



(91) 



and focus on the region very close to the shell atv = 0. Consider a point Pin (Pout) just inside (outside) the shell, and let 
the coordinate difference between Pin and Pout be AX = (Aw, Ar, Ax). Take another point M on the shell v — 0, and 
let the coordinate difference between Pin (Pout) and M be AXin — (Auin, Arin, Axi^) {AXout = (Awout: Aroutj Axout)): 
so that AXin + AXout — AX. Then the distance from Pin to Pout via M is 



A,s - ^-flAvl + 2ArinA«in + r^Ax^ + ^-f^^.Avl^, + 2(Ar - Arin)Az;out + r2(Ax - Axin)^. 
We want to find the point M that minimizes As. Extremizing this with respect to Ann, Axin, we find 



Ann Ar+\{fl~fl^,)Av,^, 



Aa 



At 



Av 



and therefore 



Arout ^ Ar+l(/,^nt-/fn)A^i, 
Az;out Av 

From these we obtain the following refraction law as required: 



An, 



AXout 
AWout 



Ax 
Av 



Ax 
Av 



dr 




dr 




^2 
' H 


dx 




dx 




dv 


in 


dv 


out 


2 ' 


dv 


in 


dv 


out 



(92) 



(93) 



(94) 



(95) 



By plugging the inside solution, (89), and the outside solution, (44) -(46), into the refraction condition (95), we 



obtain the relation between the parameters of the geodesic: 



E = ± 



27*^1, 



and J = =F- 



(96) 



where is the value of r at which the geodesic intersects the shell, and how to determine it will be explained below. 
The first sign combination applies to branch 1 for r<,^ < rn/V^ and branch 2 for > ^h/V^, while the second 
sign applies to branch 1 for rgw > rn/V^ and branch 2 for r^^ < rn/V^. Recall that branches 1 and 2 were defined 
in section 



11 A4 below eq. (46) 



It is possible that the geodesic crosses the horizon (r = r^) before reaching the shell (v — 0). When this happens, 
let us require that the geodesic goes to the boundary r = cx) at A = oo and, as we decrease A, it enters the horizon 
with v staying finite. In the previous section, we saw that the geodesic crosses the horizon at A = Xh± and, below 
(42), we observed that v is finite across A = Xh- but not across A = Xh+- Therefore, we need that A^+ < X^- 



which means that P < 0. So, we choose the second sign combination in (96) 



E= 



2r2 

sw 



and J — 



rH 



(97) 



This leaves branch 1 for Vsw > fn/V^ and branch 2 for Vsw < fn/V^- 



The equation (97) depends on r^^ which is where the geodesic hits the shell. Because the shell is at w = 0, the 
value of is determined by 



= v{rs^) 



1 

2^ 



H 



log 



(1 + E)rl ± -{l + J^- E^)rj,rl^ + J^rj, 



TH rl^ - (1 - E)t\ ± y^t; 



(1 + J2 - E'^Ym 



Jh 



H 



+ to 



(98) 
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where in the second equahty we used (46) and the ± signs correspond to branch 1 and 2, respectively. If we rewrite 
E, J in favor of r^wi f* using the relation (97 1, this becomes 

= J— log 



2rH 



(99) 



Note that the two possible signs in (98 1 lead to the same condition ( |99[ ). To solve this, let us make the convenient 
definitions 



P ■ 



p sin 9 . 



The range of the parameters is ^? € [0, p > 0, and a > 1. In terms of these, the relation (99) becomes 



4p(l + c) 



1-i 
a 



c + 2p(l + p)(l + c) 
with c = cos9. This can be solved for p, giving two branches of solutions: 



p{a,c)- 



a + 1 1 
2(a-l) 2 




(100) 



(101) 



(102) 



We can see that only p+ is allowed, as follows. First, we can show that, for fixed c, the derivative of p{a,c)± 
with respect to a is always negative (positive) in the range a £ (l,oo), c e (0,1] [69:, so that p(a, c)± monotonically 
decreases (increases) as we increase a. Furthermore 



p{a,c)± 

Therefore, for any finite values of t^, 

p{a,c)- < 




1 ± 



1 - 



as a — 7- oo (or to — ?• oo) , 



< P(a,c)+. 




(103) 



(104) 



For > riil^fi (i.e. p > 1/a/2), this means that we should take the p+ branch. The case with < rn/V^ (i-C. 
p < l/\/2) needs some more work. We only need to consider geodesies that do not reach r — 0, since otherwise we 
cannot glue branch 1 to branch 2 at Xmin before the geodesic reaches the singularity. Thus, when < rn/V^, we 
need A± < and B± > 0, that is J < 1 + i?, for < and < J < 1. Using (97), the last condition translates into 
the following condition for p: 



(105) 




This means that we should take p+ in (102) for all c G [0, 1]. 

In Figures 10 and 11 geodesies in the {v, r) plane are plotted for different values of tg and £. The geodesies start 
from the boundary at(r, v) = (oo,to), plunge in the bulk in the planar BTZ geometry, eventually refract at the shell 
at w = and propagate in pure AdS. All profiles are symmetric under x — > —x. The curve in black dashed is the 
apparent horizon. 

In Figures [12] and |13[ we display the boundary time and spatial separation dependence of the geodesies in the 
Carter-Penrosc diagram of Vaidya spacetime. The transformations to the {U, V) coordinates are given by 



U 



tan 



2 



t + 



1 



V 



tan 



for V <0, 



U = — tan 

TT 



tanh 



2 



1 ,-ir 
t H coth ^ — 



V 



tan 



2 



tanh 



2 



1 , 1 ^ 

t coth ^ — 



(106) 



(107) 
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FIG. 10: Geodesies in the (i;, r) plane for fixed to = 2 
and £ ~ 3.0 (brown, top), £ ~ 4.6 (orange, middle), 
£ ~ 68.2 (red, bottom). In black dashed, the apparent 
horizon. 
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FIG. 11: Geodesies in the {v,r) plane for fixed £ « 
21.3 and, from left to right, to = 0.1 (brown), to — 1 
(orange), to — 4 (red), to ~ 10.6 (blue). In black 
dashed, the apparent horizon. 



for V > and r > rn, and 



U = — tan 

TT 



COtll 



t H tanh ^ — 



V 



tan 



tanh 



t tanh ^ — 



(108) 



for w > and r < tr- In these coordinates, the "singularity" is located dX U + V = 1 and is depicted by a wavy line 
in Fig. 12 and 13 the shell is at F = and the horizon is at J7 = 1/2 (dashed line). 





FIG. 12: Geodesies in the spacetime Garter-Penrose 
diagram for fixed to = 2 and £ ~ 3.0 (brown, bottom), 
I ~ 4.6 (orange, middle), £ ~ 68.2 (red, top). Surfaces 
of constant r are plotted in gray. 



FIG. 13: Geodesies in the spacetime Garter-Penrose 
diagram for fixed £ « 21.3 and to — 0.1 (brown, bot- 
tom), to = 1 (orange, middle), to = 4 (red, top). 
Surfaces of constant r are plotted in gray. We do not 
display in the diagram the geodesic with to « 10.6 (in 
blue in Fig. Illl) because it would cover the horizon. 



With all the relations above in hand, we can compute the geodesic length and spatial boundary separation in terms 
of the parameters of the geodesic. The geodesic length is the sum of the geodesic length for the inside part, which 
can be computed from (53), and the one for the outside part, which can be computed from (43). The result is 



(5/: = 2 In 



-In 



4^4 



(109) 



Here, just as in (52), 5£ has been renormalized by subtracting an IR divergent quantity 21n(2ro). The spatial 
boundary separation can also be computed by summing the contribution from the inside part ( 89 ) and the outside 
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part (45), the result being 



In 



2r (r^ 
^' sw \' s 



(110) 



Although dC and £ are written in te rms of r^u, and r.^ in the above expressions, we can write them in terms of to 



and s = Vl — using the definitions ( 100 1 and plugging in the solution p — p{a, c)+ of Eq. ( 102 ) 



2p = coth(rH to) + \ coth (rnto) - 



2c 



Explicitly, 



6£{to,£) = 21n 



sinh(rHto) 



where s{£,to) € [0, 1] is parametrically defined by 



2c 
sp 



In 



2(l + c)p2 + 2s/9-c 
2(1 + c)p2 ~2sp~c 



(111) 



(112) 



(113) 



^. d > 2-dimensional field theories: numerical analysis 



As in the AdSa case, we consider a four- or five-dimensional dynamical metric which interpolates between pure 
planar AdS at early times and a Schwarzschild black brane at late times. The transition is induced by an infalling 
shell of null dust. The metric in {d + 1) dimensions was given in Eq. ([!]). We again consider the mass function Q, 
where the profile parameter is set to vq = 0.01 in the remainder. We consider geodesies with a boundary separation 
along xi, denoted x in the following, while the {d — 2) coordinates (x2, . . . ,Xd-2) of both endpoints are the same. 
This implies that z and v only depend on x, i.e. z = z{x) and v = v{x). We assume that x runs between —£/2 and 
+£/2. Then the length of the geodesic is given by: 



C = 



e/2 



2z'v' 



z{x) 



(114) 



where ' = d/dx. We notice that the integrand has no explicit x-dependence, implying the existence of a conserved 
quantity, similar to the black brane case. The conservation equation reads: 



1 - (1 - m{v)z'^)v^ - : 
The two equations of motion following from ( 114[ ) are 

zv" + 2z'v' -l + V 
m(v) 



d-2 



d„.l2 



z" + (1 - m{v)z'^)v" 



z^v'^ 



2 

d'm{v)z 



m{v)z V' 



"-^z'v' 



= 0, 
= 0, 



(115) 

(116) 
(117) 



where m{v) = dm{v)/dv. It is possible to show that (1151 and (116) imply (117) after taking the derivative of ( |115 ) 
with respect to x. Therefore we restrict ourselves to solving (115) and (116). We construct solutions z{x) and v{x) 
with the X — >■ —x symmetry which is already present in the equations of motion, and which satisfy the boundary 
conditions 



z(0)=z*, u(0) = w*, w'(0) = z'(0). 



(118) 



After constructing a geodesic for a specific choice of values (z*,w*), information about the boundary separation and 
boundary time at which the geodesic is inserted follows from 



z{£/2) = ZQ, v{£/2) = to. 



The on-shell length is obtained from Eq. (114), upon use of the conservation equation (115), 



e/2 



dx 



z(xf 



(119) 



(120) 



This length should be regularized by subtracting the cut-off dependent part 21n(2/zo), yielding 5£(^, to)- 
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3. Thermalization 



Having collected the geodesic lengths in the infalling shell backgrounds in AdS3^4_5, we can use these to describe 
the process of thermalization following a quench. To this end, we measure the approach to thermal equilibrium by 
comparing SC at any given time with the late time result (5>Cthormai- In any dimension, this compares the logarithm of 
the two-point correlator at different spatial scales with the logarithm of the thermal correlator. However, we find it 
more revealing to consider C = 6C/£, where we divide by the spatial separation on the boundary (we will discuss the 
significance of this later in Section HIC where we deal with entanglement and Kolmogorov-Sinai' entropies). In figure 
[T4]we plot this measure for two-, three- and four-dimensional field theories. In all three cases we observe a delay in 
the onset of thermalization. The reason for this "delay" is simply that the effect of the medium only becomes fully 
apparent at distances of the order of the thermal screening length ^ Although a very small volume of 

linear dimension ^ ^ would appear fully thermalized after a time to = ^/2, its contribution to the entropy of a 
large volume would be disproportionately small, because it does not support the momentum modes that constitute 
the thermal medium at large. As a consequence, the rapid linear increase of the logarithm of the two-point function 
seen in Fig. [14] thus only sets in after some delay. Obviously, this effect is more pronounced in higher dimensions. 
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FIG. 14: SC — (5-Cthormai {£■ = jC / £) sls SL functloii of boundaiy time to for d = 2 (left), d = 3 (middle), and d = 4 (right) for a 
thin shell {vo — 0.01). The boundary separations were taken to be ^ = 1,2,3,4 (from top to bottom curve). All quantities are 
given in units of M. These numerical results coincide with our analytical expressions in AdSa (d = 2) in the limit vq — >■ 0. 



From the curves in Fig. 14 we can extract different thermalization times for any spatial scale: 



(1) The critical time Tcrit at which the tip of the geodesic grazes the middle of the shell at z; = 0. This can be 
computed by asking when a geodesic with a given boundary separation in the black brane geometry outside the 



infalling shell just grazes the latter: Tcrit{£) — J. 



dz 



where is determined by the boundary separation 



(2) The half-thermalization time Ti/2- time that measures the duration for the curves to reach half of their equilib- 
rium value. 



time at which thermalization proceeds most rapidly, thus for which the curves in Fig. 14 are 



(3) The time t-^ 
steepest. 

For two, three and four dimensional field theories all these times are plotted in Fig. [15) In c? = 2 we can analytically 
derive the linear relation Tcrit = ^/2, as also observed in [?7] . 

The linearity of Tcrit {(■) in two dimensions is expected from general arguments in conformal field theory |55j , and the 
coefficient is as small as possible under the constraints of causality following a quantum quench. The thermalization 
time scales ti/2 and Tmax for three- and four-dimensional field theories (Fig. 15 middle and right) are sublinear 
in the spatial scale. In the range we study, the complete thermalization time Tcrit deviates slightly from linearity, 
and is somewhat shorter than £/2. These observations pose the question whether a rigorous causality bound for 
thermalization processes exists or not. We will come back to this question below, after examining other probes than 
two-point functions. One of the issues here is that the initial conditions are homogeneous - thus different domains will 
be independently coming to equilibrium at the same temperature possibly leading to apparent violations of causality. 

General arguments for a sharp quench in a 2-dimensional field theory [271 155j predict a nonanalytic feature where 
thermalization at a spatial scale £ is completed abruptly at Tcrit (^)- This is evident in Fig. 14 (left) as a sudden change 
in the slope at Tcrit, smoothed out only by the small non-zero thickness of the shell, or equivalently, by the intrinsic 
duration of the quench. We find a similar (higher-order) non-analyticity for d = 3, 4 (Fig. 14 middle and right) and 
expect this to be a general feature for quantum quenches in all dimensions. 
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FIG. 15: Thermalization times (rcrit, top line; Tmax, middle line; ri/2, bottom line) as a function of spatial scale for d = 2 (left), 
d = 3 (middle) and d — 4 (right) for a thin shell {vo = 0.01). All thermalization time scales are linear in i in two dimensions, 
and deviate from linearity in £ in three and four dimensions. 



Fig. 15 shows that complete thermalization of the equal-time correlator is first observed at short length scales, or 



large momentum scales (see also [25^), in contrast with the behavior of weakly coupled gauge theories. In the "bottom- 
up" scenario [S] applicable to that case, hard quanta of the gauge field do not equilibrate directly by randomizing 
their momenta, but gradually degrade their energy by radiating soft quanta, which fill up the thermal phase space and 
equilibrate by collisions among themselves. This bottom-up scenario is linked to the infrared divergence of the splitting 
functions of gauge bosons and fermions in perturbative gauge theory. It contrasts with the "democratic" splitting 
properties of excitations in strongly coupled SYM theory that favor an approximately equal sharing of energy and 
momentum |46| . One might, therefore, have expected that thermalization proceeds fundamentally differently ("top- 
down") in strongly coupled gauge theories, and this is evident within the AdS/CFT paradigm with the natural initial 
conditions used here. 



B. Equilibration of Wilson loops 

Another nonlocal probe of thermalization is the Wilson loop expectation value. As we discussed in Sec. |1IB[ this 
quantity is related in the AdS/CFT correspondence to the area of a minimal surface in AdS space that bounds the 
desired loop on the AdS boundary. In this section we compute such minimal areas in the infalling shell background. 



1. Circular Wilson loops 



We first consider minimal surfaces corresponding to circular Wilson loops in the Vaidya background. The corre- 
sponding minimal surfaces were analyzed for the vacuum (pure AdS) and finite temperature equilibrium (AdS black 



brane) situations in Sect. II B 1 This set-up was also analyzed in p8] in the context of studying entanglement entropy. 
The parameterization of the surface and boundary conditions are as given in Section |I1B 1[ the only difference being 
that the mass function m{v) is now given by (|4|. The area functional is a slight modification of ( [66| {A = a'^Nc) 
given by 



A{tQ,R)= f dp-^ y'l - (1 - m,{v)zd) u'2 _ 2z'v' , (121) 



where ' = d/dp. The explicit p dependence again means there is no conservation equation. The equations of motion 
become quite involved and we omit them here. 

We solve the equations of motion numerically, and although we want our output as a function of the boundary 
radius R, for practical purposes we have to input boundary conditions at the tip of the surface at p = where z = 
and V — v^. For each value of we consider, we find numerical solutions for the functions z{p) and v{p) for various 
values of w*, each of which reaches the boundary (which is cut-off at zq) at poiv^,). We then make an interpolation 
of po{vt), which allows us to find the exact value of such that po{v^^,) = R. For each value of z*, we now know the 
value of v^, giving the surface with boundary radius R, and we calculate the area (using the functional above) and 
the boundary time to = v{R). These are the quantities we finally plot. There are further numerical challenges with a 
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divergence as p — )■ and with integrating our solutions near the AdS boundary, where contributions to the area are 
weighted much higher. Thus we used an expansion around p = 0, and different techniques to deal with a separate 
integral near the boundary, as necessary. 

As in ( 72 ) we regulate the boundary by subtracting the cut-off dependent piece of the AdS area, 



5A{to,R) =A(to,R) 



R 



(122) 



We define SA = 6A/{-kR?) dividing by the area of the region on the boundary bounded by the loop. We then 
plot 5 A — i5.4thormai as a function of the boundary time, subtracting the black brane value to show the approach to 
thermalization. This is shown in Fig. [16] for AdS4 and AdSs for a series of values of R. 
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FIG. 16: SA — SAthcimiii {A = A/{/ttR^)) as a function of to for circular Wilson loop radii R — 0.5, 1, 1.5, 2 (top curve to bottom 
curve) and mass shell parameters vo ~ 0.01, M = 1, in three-dimensional (left panel) and four-dimensional (right panel) field 
theories. 



Repeating the analysis performed for Wightman functions, we calculate the three thermalization times defined 
earlier (see Sec. Ill A3) as a function of the loop diameter (Fig. 17 1. We use the diameter rather than the radius here 



as it is the analogue of the separation £ that we plotted for the geodesies. The complete thermalization time Tciit(-D) 
is close to being a straight line of slope 1/2 for three dimensional theories over the range of scales that we study 
(also see PS^) (it would be unit slope as a function of the radius). But in four dimensions Tciit{D) deviates somewhat 
from linearity and is shorter than D/2. Overall, our thermalization times for Wilson loop averages are remarkably 
similar to those for two-point correlators . This suggests that higher-order correlators thermalize similarly to the 
basic Green function in these strongly coupled theories, so that any of these non-local probes fairly assesses progress 
towards thermalization. 



1. 




1. 




FIG. 17: Wilson loop thermalization times (rcrit, top line; Tmax, middle line; ri/2, bottom line) as a function of the diameter 
for circular Wilson loop operators in three-dimensional (left) and four-dimensional (right) field theories. 
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2. Infinite Rectangular Strips 



As a second example we consider an infinite strip, similarly to (28j. A key difference between the strip and the 
Wilson circle is that the size of the latter is set by a single scale, the radius. The strip could be regarded as a limit of 
a ellipsoidal loop with a highly elongated semi-major axis. Thus there are in effect two scales involved in the Wilson 
strip - the width and the (infinite) length. As before, at early times the background is AdS4, but evolves to a black 
brane at late times due to an infalling shell of null dust. The metric is at all times given by ([I]) with c? = 3 or 4. 
Again, the mass function m(v) determining the evolution of the metric as a function of the bulk light-cone time v is 
modeled by the expression Q. We set the parameter vq = 0.01. At the AdS4 boundary we consider a rectangular 
strip parametrized by the coordinates xi and X2, such that 



XI e (-^/2,^/2), x2eiO,R), 



(123) 



where R will be taken to infinity [70]. The coordinate xi is denoted x in the rest of the section. The area of the string 
surface with this rectangular base is given by 



Aito,£,R) 



R 
2^ 



e/2 



(1 - m{v)z'i)v''^ - 2z'v' 



when we consider the following embedding profile for the string: 

v = v{x), z = z{x), z(±^/2) = zo, i;(±£/2) = to. 



(124) 



(125) 



We notice that z and v only depend on x and that x itself does not appear in the action, similarly to the geodesic 
length in a dynamical Vaidya background. There exists a conserved quantity corresponding to the maximum value 
of z denoted as z*. The conservation equation then reads: 



1 - (1 - m{v)z'^) 



2z'v' 



From the action ( 124 ) we obtain the following two equations of motion: 

d-4 

m{v) 



+ Az'v'-2 + 2v''^ + m{v)z'^v''^ = 0, 



z" + (1 - m{v)z'^)v" 



-z^v'^ 



dm{v)z'^'^z'v' = 0. 



(126) 

(127) 
(128) 



One can show that taking the de rivat ive of ( 126[ ) with respect to x and combinig with (127) leads to (128 1. This 
it is sufficient to solve (126) and ( |127[ ). In order to do this numerically we use the symmetry of the string surface, 
i.e. we construct one half of the solution z{x) and v{x) and then use the reflection symmetry to construct the other 
half. Our boundary conditions are 2(0) = z*, v{Q) = u*, v'{Q) = = z'(0). After constructing a string surface 
satisfying these boundary conditions, we read off the boundary separation and the boundary time using z{£/2) = zq, 

v{e/2) = to. 



Using the symmetries of the embedded string surface and the conservation equation (126), we obtain the on-shell 
area as: 



R 



e/2 



A{to,e,R) = - / dx- 

TT Jo ' 

We regulate the area by subtracting the cut-off dependent part 

6AitoJ,R)^A{to,£,R) 



1 R 
Zf) vr ■ 



(129) 



(130) 



We then compare the area in the shell background with the area in a thermalized background at different times for 
flxed scales in Fig. 
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where we defined SA = SA/{Ri/n). In four dimensions {d = 3) it turns out that there is a 
range of times to and Wilson strip widths £ for which there are three different minimal surfaces in AdS4. Thus, in this 
case, we observe a "swallow tail" in the Wilson loop thermalization curve for large enough boundary separations, as 
was also noted in I28| . However, in five dimensions (d = 4) the thermalization curves do not exhibit the swallow tail. 

Ti/2 and T„ 



In Fig. 
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we plot the three different thermalization times r^^tj T1/2 and Tmax for d — 3 and d — 4 dimensional field 
theories, in three and four dimensions, the time for complete thermalization, Tcrit (dashed green curve), exceeds the 
linear relation £/2. The half-time T1/2 (orange curve) is sub-linear. The time r^ax (dot-dashed blue curve) exhibits a 
different behavior in three and four dimensions. The time for complete thermalization for a strip of width £ is slower 
than the thermalization time of a circular loop of radius R. This may be because the long direction of the strip 
provides a second, larger scale. 
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FIG. 18: SA — 5.4thcrmai (A = A/{R£/it)) as a function of boundary time to for d = 3 (left) and d = 4 (right) for a thin shell 
{vq = 0.01). The boundary separations were taken to be ^ = 1,2,3,4 (from top to bottom curve). All quantities are given in 
units of M. The middle panel shows a zoomed-in version of the swallow tail for £ — i, where we plot 1{SA — 5.4thcrmai) to 
amplify the effect. 




FIG. 19: Wilson strip thermalization times (rcrit, top line; Tmax, middle line; ri/2, bottom line) as a function of spatial scale for 
d — 3 (left) and d = 4 (right) for a thin shell (vq — 0.01). All thermalization times deviate from linearity in the spatial scale. 



C. Equilibration of the entanglement entropy 

Finally we consider how the entanglement entropy equihbrates following a quench. Because it measures all contri- 
butions to the information loss caused by the restriction of the field theory to a finite volume, this quantity provides 
for a more comprehensive measure of equilibration than either the two-point function or the Wilson loop expectation 
value. The entanglement entropy can be considered as a special case of the standard coarse grained entropy of equili- 
brating quantum systems |58j . As discussed in Sec. |IIC the entanglement entropy in a two-dimensional field theory is 



related to geodesic lengths in AdSa, while in three-dimensional field theory entanglement is related to minimal surface 
areas in AdS4. Both of these quantities have been computed above in the infalling shell background in the process of 
studying two-point functions and Wilson loops. It is interesting that the Wightman function in 2 dimensions and the 
Wilson loop in 3 dimensions are so closely related to entanglement. 

For 4-dimensional theories, to study the entanglement entropy in spherical regions, we need to compute the volume 
of minimal three-surfaces with S'^ cross-section in the asymptotically AdSs infalling shell geometry. As with the pure 



AdS and black brane cases we discussed in Section II C the method is a straightforward generalization of the circular 



Wilson loop case, so we do not repeat the details of Section pi B 1[ We replace the black brane tension M in equation 
(82 1 with the dynamical mass function (Eq. (|4|). We can regulate the resulting volumes by subtracting the divergent 



part of the 3- volume of the 3-surface of the same boundary radius in pure AdS, which is given analytically in (83). We 
call this regulated volume SV. However, as in previous cases, to illustrate the thermalization process we find it more 
instructive to subtract the renormalized 3-volume of the solution of the same radius in the black brane background. 



SVthermai ■ We ultimately plot SV — SVthermai against the boundary time in the left panel of Fig. 20 V is defined as 
V divided by the volume of the region that it bounds on the boundary. As in previous cases we can also calculate the 
three thermalization times we defined earlier for the 3-volume, at different values of the diameter D = 2R. These are 
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shown in the right panel of Fig. ( 20 1 . 

0. 




b- 1. 




to 

FIG. 20: The left panel shows 5V — <5Vthormai (V = V/(47rii^/3)) as a function of boundary time tg for a 3-surface in d = 4 for 
a thin shell {vo — 0.01) and 4 equally spaced boundary radii R — 0.5, 1, 1.5, 2. The right figure shows the thermalization times 
{Tcrit, top line; Tmax, middle line; ti/2, bottom line) as a function of spatial scale for the same 3-surface. 



Our thermalization times for Wilson loop averages and entanglement entropy seem remarkably similar to those 
for two-point correlators. Slightly "faster-than-causal" thermalization, possibly due to the homogeneity of the initial 
configuration, seems to occur for the probes that do not correspond to entanglement entropy in each dimension. For 
the latter, the thermalization time is linear in the spatial scale and saturates the causality bound. As the actual 
thermalization rate of a system is set by the slowest observable, our results suggest that in strongly coupled theories 
with a gravity dual, thermalization occurs "as fast as possible" at each scale, subject to the constraint of causality. 

The average growth rate of the coarse grained entropy in nonlinear dynamical systems is measured by the 
Kolmogorov-Sinai (KS) entropy rate /iks [591, which is given by the sum of all positive Lyapunov exponents. For a 
classical SU(2) lattice gauge theory in 4 dimensions, /iks has been shown to be proportional to the volume [60J. For a 
system starting far from equilibrium, the KS entropy rate generally describes the rate of growth of the coarse grained 
entropy during a period of linear growth after an initial dephasing period and before the close approach to equilibrium 
[61j . Here we observe similar linear growth of entanglement entropy density in d = 2, 3, 4 (leftmost panels Figs. 14 16 
|20[). For small boundary volumes, the growth rate of entropy density is nearly independent of the boundary volume 
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(almost parallel slopes in the leftmost panels of Figs. [I4j [l6j [20| and nearly constant maximal growth rate in Fig 
left). Equivalently, the growth rate of the entropy is proportional to the volume - suggesting that entropy growth 
is a local phenomenon. However, in d = 2 where our analytic results enable study of large boundary volumes £, 
we find that the growth rate of the entanglement entropy density changes for large falling asymptotically as X/l 
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b panel), and thus cannot arise from a local phenomenon. This behavior suggests that entanglement 
entropy and coarse grained entropy have different dynamical properties. 

In summary, we have investigated the scale dependence of thermalization following a quench in 2-, 3-, and 4- 
dimensional strongly coupled field theories with gravity duals. We found that the entanglement entropy sets a time 
scale for equilibration that saturates a causality bound. Our results raise interesting questions about the relationship 
between the entanglement entropy growth rate and the KS entropy growth rate defined by coarse graining of the 
phase space distribution. 



IV. DISCUSSION 



In this work, we explored the approach to thermal equilibrium of simple nonlocal observables in strongly coupled 
conformal field theories with a gravity dual. In particular, we investigated the behavior during thermalization of 
equal-time two-point functions, Wilson loops and entanglement entropy in boundary field theories holographically 
represented by asymptotically AdS geometries. The thermalization process was modeled in the bulk by the collapse 
of a spatially homogeneous thin shell of null dust that eventually forms a black hole. The thickness vq of the shell is 
related to the duration of the process of energy deposition in the boundary field theory. In the limit uq <C , where 

parametrizes the effective UV cut-off of the field theory, the energy deposition can be considered as instantaneous. 

Our observables were calculated in a semi-classical approximation appropriate for high-dimension operators, in 
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/ or D / / 

FIG. 21: (Left) Maximal growth rate of entanglement entropy density vs. diameter of entangled region for d = 2,3,4 (top to 
bottom). (Middle) Same plot for d — 2, larger range of £. (Right) Maximal entropy growth rate for d = 2. 

which they correspond to geodesies or minimal surfaces in the bulk. For the AdSs case, where the thermal limit is 
represented by the planar BTZ black hole, we were able to solve the problem analytically. Not only did this provide 
a check on our numerical approach, which is required for higher dimensional cases, it also allowed us to obtain results 
in regimes beyond reach of our numerical analysis. 

In all cases we found that the thermal limit is reached after a finite time Tcrit, which is a function of the geometric 
size of the probe in the boundary field theory, e.g. the separation £ of the two points of the equal-time Wightman 
function or the radius R of the circular Wilson loop. For those cases where the logarithm of the correlation function 
is proportional to the entanglement entropy of the enclosed area on the boundary, i.e. for geodesies in AdSa, minimal 
surfaces bounded by a circle in AdS4, and minimal co-dimension two surfaces bounded by a spherical shell in AdSs, 
we found that Tcrit — (AdSa) and Td-a = D/2 (AdS4,5), in the regimes we were able to study. Our result confirms a 
general rule for two-dimensional conformal field theories implied by causality (as also observed in |27j ) and generalizes 
it to higher dimensions (as also done in for AdS4). 

The conclusion we can draw from these results is that the decoherence and equilibration of instantaneously deposited 
energy generally propagates at the speed of light in the conformal field theory. This makes sense, because the process 
can be qualitatively understood as a cascade of gauge boson splittings, which transfers energy from all momentum 
scales into the infrared until all energy is distributed thermally. The speed by which this process evolves is determined 
by the speed of propagation of the gauge quanta, i.e. the speed of light. Alternatively, if the initial energy deposition 
were spatially inhomogeneous, one would expect the dissipation of the local inhomogeneities in the initial energy 
density to be constrained by the speed of sound, which is equal to c/ -^/S in a conformal field theory. This argument 
invites the conjecture that local equilibration of energy, in the sense of hydrodynamics, is controlled by the speed of 
light in conformal field theories, whereas global equilibration requiring hydrodynamic transport of energy is governed 
by the speed of sound. On the other hand, we found that the naive causality argument does not hold for the two-point 
function in 3d and 4d field theories and for Wilson loops in 4d field theories, which do not have an interpretation as 
entanglement entropy. In this case we obtained Tcrit < ^/2. This is not as surprising as it might sound: we already 
know of other observablcs that thermalize faster than the entanglement entropy, for instance one-point functions. 

We also found that the transition to full thermal equilibrium is abrupt and non-analytic, as expected from the 
causality argument. One might have objected that the non-analyticity of the transition to full thermal equilibrium is 
probably due to the geodesic approximation and will be smoothed out once we consider the exact Wightman function. 
Corrections to the geodesic approximation can be worked out by considering the exact first-quantized path integral 
representation of the propagator, whose saddle points are given by geodesies. Possible corrections are then due to 
fluctuations around the saddle point. However, if the geodesic saddle-point lies entirely within the AdS Schwarzschild 
geometry and does not graze the shell, such higher order corrections will not be sensitive to the presence of the shell, 
since the shell is invisible in perturbation theory around the geodesic. This remains valid when we couple the heavy 
bulk scalar field to other bulk degrees of freedom. At best, there could possibly be non-perturbative corrections due 
to other complex saddle points which are sensitive to the existence of the shell, but is unclear how to detect those in 
time-dependent situations. Therefore, for sufficiently small separations the two-point function appears to be exactly 
given by the thermal answer, which will then abruptly change to a non-thermal answer once the geodesic starts 
crossing the shell. As a result, we expect the non-analyticity to be valid beyond the geodesic approximation. The 
non-analyticity is probably a consequence of our treatment of the endpoints of the Wightman function, which are 
both located at the same bulk variable zq. The origin of a similar problem has been studied in the theory of boundary 
critical phenomena, where one typically introduces an extrapolation length tq to tame the sharp UV cut-off of the 
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conformal field theory [551 This corresponds to the introduction of a smooth UV cut-ofF of the form ex.p{—ETo), 
where E is the energy. The sharp cusp in the asymptotic result is then rounded over a region |< — £/2| ~ tq. In our 
holographic setting, a similar approach would smear the two endpoints of the correlation function independently over 
some range of values zq. We have not studied the consequences of such a prescription. 

One may wonder whether our results are in tension with those of I64j . where instantaneous thermalization was 
found for two-point functions related to Brownian motion of a "quark" represented by a string stretching from the 
boundary to the horizon of an AdSs-Vaidya spacetime. While the two-point functions studied in |64j did not involve 
high-dimension operators, so that there is no reason to expect a geodesic approximation to be valid, one can get 
some intuition by computing geodesies in the induced geometry on the string worldsheet. We have done so, and 
find that geodesies connecting two points on the string outside the shell never cross the shell. So within a geodesic 
approximation, this would explain why the observable considered in |64| thermalizes instantaneously, in contrast to 
the observables studied in the present paper. 

The fact that correlation functions of a small geometric size thermalize earlier than those of a large size implies 
that high momentum modes in the boundary field theory approach thermal equilibrium faster than long wavelength 
modes. In other words, the boundary field theory thermalizes first in its ultraviolet domain and later in the infrared 
or, as one might say, thermalization proceeds from the top down. In part, this is a consequence of the fact that 
the specific energy injection mechanism considered here has support in the UV: the energy shell is injected into the 
geometry near the boundary and then falls into the deeper regions of the AdS space, which represent IR modes of the 
boundary field theory. While this appears natural when viewed from the gravity side, it represents a radical deviation 
of the thermalization behavior of the dual field theory from that known in weakly coupled non-Abelian gauge theories, 
where thermalization occurs from the bottom up, independent of how the energy is injected into the field theory. For 
the reader unfamiliar with heavy ion collision phenomenology, we point out that the scale at which energy is deposited 
in the collision - the parton saturation scale Qs of the colliding nuclei - is higher by an inverse power of the coupling 
constant as than the thermal scale T after thermalization. Thus, while the energy is not injected in the extreme UV 
in a nuclear collision, it is certainly not injected in the infrared. We also note that the difference in thermalization 
behavior is closely related to the observation that highly localized excitations in the strongly coupled gauge theory 
do not evolve into jets as is commonly found in weakly coupled, asymptotically free gauge theories that permit a 
perturbative treatment. 

It is tempting to speculate to which extent our results may apply to the thermalization of QCD matter that is 
produced in relativistic heavy ion collisions. Phenomenologically, it is known that the thermalization has to occur 
fast, allowing the matter to expand according to the laws of nearly ideal hydrodynamics. If the relevant length scale 
for thermalization is given by the thermal scale £ ~ Ji/T, our results suggest that Tdit ~ 0.5h/T if the matter is 
strongly coupled. For initial temperature value T 300 — 400 MeV at heavy ion collider energies, we obtain the 
estimate Tc^it ~ 0.3 fm/c, comfortably short enough to account for the experimental observations. We know that such 
a short thermalization time applies to all momentum scales, because jets are created in the nuclear collisions, albeit at 
suppressed levels. However, as already mentioned, the relevant momentum scale for energy deposition is the parton 
saturation scale Qs of the colliding nuclei, i.e. the transverse momentum below which most partons in the colliding 
nuclei are found. For achievable energies and large nuclei, Qs < 3 GeV; in this range of momenta QCD may well be 
considered as a strongly coupled gauge theory. 

Though this was not the emphasis in this paper, the rapid thermalization studied in this paper may also be 
of relevance for black hole physics. It has been suggested that black hole creation is the fastest possible form of 
thermalization that exists in nature [6 5) and it would therefore be of particular interest to more directly link the 
causality bound to this conjecture. 

We also encountered phenomena that do not appear to have a simple explanation. For example, for long rectangular 
Wilson loops in the boundary gauge theory of AdS4, the area of the minimal surface exhibits a "swallow tail" behavior: 
the transition to the thermal limit proceeds via a discontinuous succession of shapes of the minimal surface. It is 
unclear why this phenomenon occurs for the rectangular Wilson loop but not for the two-point function, as the 
differential equation for the stationary surface is quite similar to that for the geodesic curve. 

Finally, we briefly comment on the difference between dynamically collapsing shells and quasistatic shells that 
adiabatically approach the event horizon. As we show in the Appendix, the swallow tail behavior is ubiquitous for 
quasistatic shells close to the horizon, leading to significant differences between dynamical and quasistatic shells for 
large boundary separations. 
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Appendix A: The quasistatic approximation 

In the main text we considered a dynamical shell of matter falling into AdS space as a model of thcrmalization. An 
adiabatic approximation to the process of thcrmalization can be constructed by treating the computation of correlation 
functions and Wilson loops in a quasistatic manner. To do this we observe that the dynamical shell is falling into AdS 
space on a null trajectory. Then, at any given time, we imagine a shell located statically at the corresponding location 
on the null trajectory and compute Wightman functions, Wilson loops and entanglement entropy via the geodesic 
and minimal surface methods described in the previous section. In this quasistatic approximation the geodesies and 
minimal surfaces with endpoints at a given boundary time remain localized on the equal time surface in the bulk of 
AdS space (Fig. [2]). 

As we will see, in the quasistatic approximation, there are generically multiple geodesies connecting a pair of 
endpoints on the AdS boundary. We demonstrate this explicitly for AdSs and have found similar results for AdS4 
(not shown). There can also be multiple minimal surfaces that trace out a given loop on the AdS4 boundary. This 
"swallow tail phenomenon" was also described in the full dynamical setting for the strip Wilson loop (see [28] and 
above). When multiple geodesies or minimal surfaces are present in a dynamical setting, it will generally be necessary 
to analyze which of these saddlepoints lie within steepest descent integration contours in the path integral. Here, for 
completeness, we document the situations we have found where multiple saddlepoints exist (see below and Table |l|. 
Note that in a quasistatic setting the saddlepoint of minimal action will dominate. 



holographic probe 


quasistatic 


dynamical 


geodesies in AdSs 




X 


geodesies in AdS4 




X 


strip in AdS4 






circular loop in AdS4 




X 



TABLE I: Overview: swallow tail (-y/); no swallow tail (x) 



1. Geodesies in the quasistatic approximation 



We study the equal-time geodesic in the shell background in the quasistatic approximation for d = 2. In the 
quasi-static approximation, the AdSs geometry with a shell at a fixed radial coordinate r = Ts can be written as 



' r^dX^ 



/(' 



l/out(r)2 = r^ 



,.2 



r < rs 
r > Ts 



(Al) 



where th (< rs) is the position of the would-be horizon. Outside the shell, r > r^, the metric ( Al I is the planar black 
brane metric that we studied in the text in |II A 4[ while inside the shell, r > rs, it is t he e mpty AdSa we discussed 



around (53)-(55). Therefore, for studying the equal-time geodesic in the shell spacetime ( Al I, we can use our previous 



results. Outside the shell, the relevant expressions are found in (34), (351, (41 1, and (45 1. The outside geodesies are 



parametrized by two parameters E and J which correspond to the energy and angular momentum, respectively. We 
take i? = so that they describe equal-time geodesies (see below (48)). Also, we can take J > without loss of 
generality, because flipping the sign of J just corresponds to the reverse parametrization of the same geodesic. Inside 
the shell, the expressions for geodesies are found in ([53|, (54), and (89). The inside geodesies are parametrized by 
= 2/(!.. As is clear from (89), is the minimum value of r on the geodesic. 
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The geodesic in the entire shell background spacetime (All is obtained by gluing outside and inside geodesies across 



the shell using a refraction law so that the geodesic length is extremized. In the current situation the refraction law 
can be derived in very much the same way as it was in the Shockwave background around (95), and the result is 



fn 



dX 
dr 



fo 



dX 
dr 



(A2) 



By substituting the expression for the outside geodesic (45) and the inside one (89) into (A2), we see that the 
parameters in the inner and outer regions are related by 



r.^. = rnJ- 



(A3) 



Just as we did in the main text for the Shockwave geometry, we can compute the renormalized geodesic length 5C 
and spatial boundary separation I in terms of the parameters J, of the geodesic. The geodesic length is the sum of 
the geodesic length for the inside part, which can be computed from (53), and the length for the outside part, which 
can be computed from (43) with E = 0. The result is 



2^rl - rjjj^ 2 
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J2 



(A4) 



Here, we regularized C by subtracting the divergent quantity 21n(2ro) just as we did in the main text. Note that (A4) 
is valid if the geodesic is partly inside the shell, which is the case if J < rg/rn- In this case, (A4) gives implicitly the 
relation between i and 5C through the parameter J (note that J and are related by (A3)). 



If J > Ts/rn, instead, the entire geodesic is outside the shell and the relation between ^ and 5C is given by (52) 




6C 



6 - 




jru = 1.001 



FIG. 22: The &C, versus i curves (solid lines) for the values of rs/rn (A) sufficiently larger than 1 and (B) sufficiently close to 
1. We set r_r/ = 1. For comparison, we plotted 5C{1) for pure AdSs in dotted black lines and for AdSa Vaidya in dashed yellow 
lines. Note that most of the Vaidya curves (dashed yellow) are overlapping with quasi-static curves (solid, in various colors). 
For the values of to used for the Vaidya curves, see text (Eq. (A6l). In (B), we used (partly) different colors for the different 
regimes (i)-(iv) explained in the text in Eq. ( A5 I. In the insets are the actual shapes of the geodesies for some selected values 
of I, in different colors corresponding to different branches, in the {X,z = 1/r) coordinates (actual values of the coordinates 
not shown). The two horizontal lines in each inset represent the boundary [z = 0) and the shell (zs = l/^'s). 



In Fig. 22 we plotted 5C versus £ for different values of Ts/th- We see that the curve has more nontrivial structure 
for rg/rn ~ 1- So, let us focus on this case. If rs/rn « 1, then from (A4) one can derive that, depending on the value 
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of J, the relation between ^ and 5C has the following regimes: 
(i) ^ » 1 : ^ « ^ 



(ii) J > ^ and J - 1 

(iii) J < ^ and J - 1 



(iv) 



> J>0 







, 5C 


In 




, 5C 






In 




, 5C 




rs-TH 






, 5£ 



2 In 



2 In I + 2 In 



(A5) 



The geodesic does not cross the shell in regimes (i) and (ii) while it crosses the shell in regimes (iii) and (iv). In 
regime (i), the £-6C relation is the same as the one for empty AdSs, because the geodesic is near the boundary of 
the AdS space where the background metric is same as the empty AdS metric. In regime (ii) the geodesic length SC 
grows linearly with the boundary separation i, because the geodesic extends almost parallel to the shell just above 
it. This happens because the shell at r = is very close to the would-be horizon at r = ru- Regime (iii) is similar 
to regime (ii) except that a small middle part of the geodesic dips into the inside geometry. In regime (iv), the i-6C 
relation is the same as the one for empty AdS, up to a constant shift. This is because most of the geodesic is inside 
the shell where the metric is equal to the empty AdS metric. However, the part of geodesic outside the shell still sees 
the BTZ metric and leads to the constant shift that survives even for f — cxi. 

So, if r^/rji ss 1, within a certain range of the spatial separation there are three possible geodesies and three 
corresponding values of (5£(£); see Fig. 22 The upper bound of the range is ^max = (l/r/j) ln[(rs+r//)/(rs— r//)] while 
the lower bound can be computed from the condition dijdJ = and is given in the ru limit by Imin ~ 3.7/r^f . 
As Ts rn, we can make the upper bound Imax arbitrarily large and the "spike" in Fig. |22| arbitrarily long. 

In general, in the presence of multiple geodesies in Lorentzian signature, which geodesic makes the dominant 
contribution to the correlation function is a subtle issue |39H41]. However, in the present approximation where the 
spacetime is static, it is simply the shortest geodesic that contributes most to the correlation function. Therefore, in 
spite of the existence of the peculiar "linear" regimes (ii) and (iii), it is (i) and (iv) that determine the correlation 
function. 

For comparison, in Fig 



22 



The relevant 
at which the 



^ we also plotted SC{£) for the Vaidya geometry we studied in IIIAl 
expressions are found in (102 1 and ( |112 1. The Vaidya geodesies are parametrized by to, the time 
boundary operators are inserted. In the plot, we determined the value of Iq in terms of rg so that the shock wave in 
the Vadiya geometry is at r = at time to just as in the quasi-static metric ( Al I. Namely, because the shock wave 



is at u = 0, by setting u = in (29), we obtain the relation between to and 



2rH Ts 



Note that, in the large £ limit with fixed, the quasi-static computation (A4) gives the asymptotic behavior 



5£ w 2 In - -I- 2 In - 
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(quasi-static). 



(A6) 



(A7) 



where in the second expression we used (A6). In the same limit (which corresponds in ( 112 ) to c — > 1, r^/rn — > 
with fixed to), the Vaidya result (112) gives the following asymptotic behavior: 



(5£«21n--ln 



r7^ = 21n^+41nfcosh^') 

(V^ + i)4 2 V 2 y 



(Vaidya) , 



(A8) 



where we used a = e^^"*". These lead to the difference in the £ — >■ oo behavior of the quasi-static and Vaidya curves 
that we can observe in Fig. [22} 



2. Wilson loops in the quasistatic approximation 

In four dimensions the entanglement entropy is related to the area of the space-like string surface whose base is a 
Wilson loop at the AdS boundary. In [28J it was noticed that there can exist three different string surfaces for the 
same boundary time and the same boundary separation in a dynamical background with a thermal quench (see also 



III B 2 1 . In this subsection we will investigate whether the quasi-static approximation for Wilson loops and Wilson 



strips exhibits multiple string surfaces for the same boundary. 
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a. Circular loops 



Similar to the geodesic analysis, the four dimensional background can be divided into a black brane part (outside 
the shell) and a pure AdS part (inside the shell). However, we cannot solve the differential equations analytically and 
instead we rely on numerical methods. Outside the shell the solution satisfies the differential equation ( |68[ ) in the 
black brane background (with d — 3), while inside the shell we must find a minimal surface in pure AdS (Te. we set 
the mass M to zero in eq. (68 1. The minimal surface in pure AdS consists of a hemisphere. In the thin-shell limit, at 



the position of the shell z — z^, the minimal surface should satisfy the refraction condition 



dzo 



dp 



Vl - Mzf 



dZin 

dp 



(A9) 



This is derived analogously to the geodesic case ( A2 1 in which case it agrees with numerical results obtained by 



minimizing the path of the geodesic. Note that in the dynamical case we do not use such a condition and solve the 
equations of motion in the whole space with mass function Q. However, here, taking advantage of this simple relation 
allows to use much less computationally intensive numerics in the limit of a thin shell (where the mass function is 
a step-function). We start by constructing a hemispherical string surface whose tip is given by z{0) = z^. > Zs and 
z'(0) = 0, thus lying inside the shell. For this solution we can determine the intersection point pi with the shell and 
the gradient at that point analytically. We can then use the refraction condition to generate the rest of the solution 
in the black brane background, using the boundary condition z{pi) — Zg and z'{pi) — z'^^^. Where this part of the 
solution intersects the boundary (cut-off at zq) gives us the boundary radius R of the Wilson loop. We can also 
construct surfaces wholly outside the shell, which are just the thermal Wilson loops from section HB 1 Thus we get 
R{Zi,) for any shell position Zg. 

As with the geodesies case we can get multiple surfaces with the same boundary radius and shell position (such as 
in Fig. 22). An easy way to check whether these multiple solutions exist before calculating any areas (see also [55j) 



is to invert our numerical data appropriately so that we can plot it on a graph of against boundary time to, the 
boundary time being obtained as a function of the shell position z^: 



to 



dz 



1 - Mz^ ' 



(AlO) 



The multiple solutions will show up as multiple possible values of z^ for the same Iq. We see this in Fig. 23 ^A) where 
we have plotted for three different boundary radii. 

The area of the minimal surface consists of three parts: we add the areas of the piece in empty AdS (inside the 
shell) and of the piece in the black brane background (outside the shell) and subtract the area of the solution in empty 
AdS with the same boundary radius, 



/•Pi(to) p r 

5A{R,h)= / dp^Vl + 2'2+ / 

Jo 2 Jo 



, P 



Pi(to) 



1 



1 - Mz^ 



(All) 



In Fig. 23 ^B), we plot SA{R,to) as a function of R for a fixed shell position. The three branches appear in a very 
similar fashion to those for geodesies in AdSa shown in Fig. 22 Fig. 23 C) then shows SA{R,to) as a function of the 
boundary time to for three different boundary radii. The appearance of the three branches gives a 'swallow tail' as 
found in |28j in the dynamical case for the infinite rectangular strip, and as we also find in Section HI B 2 However, 
there are no such multiple solutions in the dynamical Wilson loop case as shown in [55| and Section HI B 1 This can 



be taken as an indication of the limits of the quasistatic approximation in accurately capturing the behavior as the 
shell gets closer to the horizon (recall that to get the three branches in Figure [23|^B) we have to wait till the shell 
approaches the would be horizon, here M^^^Zs = 1/1.001). The appearance of three branches does not, however, 
signify that something is unphysical. As mentioned, they also appear in the dynamical case for the infinite strip. 
What it does mean is that since we are using a saddle point approximation to calculate field theory quantities and our 
space time is Lorentzian, when there are multiple geodesies we should carefully follow a steepest descent procedure. 

The numerics are hard to pursue around the earliest time where the three branches occur in figure 23 [G) and we 
see some noise around one corner of the 'swallow tail'. Here the tip of the sphere at z*, which is our initial condition, 
is very close to the shell, and very small variations in give rise to large variations in R and hence SA. 



h. Strip 



We can perform a similar analysis for the Wilson strip in a quasi-static background. The equation (76) is now 
the differential equation we solve numerically in the black brane background (with d = 3). Considering a similar 
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FIG. 23: (A)z, against to for R = 1.4, 1.6, 1.8. Three ^.-values for a certain range of to makes it clear we have multiple 
branches. (B) SA against R for M^''^Zs = 1/1.001. We have multiple solutions for the same radius in the region with the spike. 
(C) SA against to for R — 1.4, 1.6, 1.8, we have three branches and the behavior of the curve there is complicated. 



boundary setup for the strip as in section [II B 2[ we follow a similar procedure to construct solutions as for the quasi- 
static Wilson loops. Also the refraction condition through the shell (in the thin shell approximation) is the same as 
for circular Wilson loops. Hence, we will not work this case out explicitly will immediately give the results. The 
string surface area consists of three different pieces: the part of the solution inside the shell (empty AdS) AaiIS, the 
part of the solution outside the shell (black brane background) Abb : f^nd a subtracting part corresponding to the full 
solution in empty AdS with the same boundary separation £. The full area in the shell background is given by 



dA{£,to) = AAds{zsito),S:) +ABB{S:,i) ~ AAds{za,e), 



(A12) 



where x represents the space-like separation of the part inside the shell. The area formula in empty AdS is given by 
the relation 



AAdsiz,y) 



R (\ 



1 ^r (-i/4)r(3/4) \ 

r(i/4)2 ; 



22/ 



(A13) 



where R corresponds to the height of the Wilson strip at the boundary (see II B 2 ) , and the black brane area is given 

by 



ABB{yuy2) = 



R 



V2 



dx- 



1 



yi 



z{x) 



z'{x)^ 
1 - Mz^ ' 



(AM) 



In Fig. 24 ^A), we plot the midpoint of z{x), namely z^,, as a function of the boundary time to- In the middle panel 
of Fig. 24 we plot the area dA{£, to) (with R/n divided out) as a function the boundary separation £ for a fixed shell 
position Zs (thus for a fixed boundary time to)- On the right panel of Fig. 24 we plot the area SA{i,to) (with R/tt 
divided out) as a function of the boundary time to for three fixed boundary separations. As in the dynamical case, 
we have three branches of solutions appearing, but here they show up already for smaller values of the radius. 



6:n 



■ m 




0.5 ^ 



(A) 



(B) 




FIG. 24: (A) 2* against to for I = 1.4, 1.6, 1.8. We see indeed the existence of three different solutions for the same boundary 
time to- (B) M/(i?/7r) against I for a fixed position of the shell: M^I'^Zs = 1/1.001. (C) M/(-R/7r) against to for i = 1.4, 1.6, 
1.8. 
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Hence, in distinction to the matter shells in early work, the shell has tension which in part contributes to the collapse 
dynamics. 

Working with renormalized operators, instead of bare operators, ensures finite expressions when the cut-off zq is removed, 
since renormalized operators do not depend on the cut-off zq. A useful heuristic treatment can be found in e.g. |50| . 



Since we are taking the coordinate x to be noncompact, in terms of the metric the spacetime (291 is identical to empty 
AdSa. Furthermore, r = Q, which is a Milne-type singularity for compact x, is not a singularity. However, we still refer to 
r = as the "singularity". As the geodesies we consider in the current paper never reach r = 0, this does not matter to us. 
dp{a, c)- /da actually vanishes when c — 0, while dp(a,c) + / da < also for c = 0. 
This is exactly the same configuration as in ,28) . 



